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We  consider  a  computer  communication  system  with  many  geographically  distributed 
stations  accessing  a  common  channel.  Stations  receive  only  channel  feedback  information. 
The  problem  is  to  design  efficient  algorithms  for  dynamically  allocating  the  channel 
resource  among  the  stations.  The  main  performance  measures  of  such  algorithms  are 
throughput  (.channel  utilization),  induced  delavs,  and  feedback  error  sensitivity. 

For  many  communication  applications  with  time  constraints  (e.g.,  transmission  of 
packetized  voice  messages),  a  critical  performance  measure  is  the  percentage  of  messages 
which  are  transmitted  within  a  given  amount  of  time  after  their  generation  at  the 
transmitting  station.  In  this  talk  we  present  a  random  access  algorithm  (.RAA)  suitable  for 
time-constrained  applications.  Performance  analysis  demonstrates  that  significant  message 
ie  lay  improvement  is  attained,  at  the  expense  of  minimal  traffic  loss. 
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We  also  consider  the  case  of  noisy  channels.  The  noise  effect  appears  at  erroneously 
observed  channel  feedback.  Error  sensitivity  analysis  shows  that  the  proposed  random 
access  algorithm  is  insensitive  to  feedback  channel  errors. 

Window  Random  Access  Algorithms  (RAAs)  are  considered  next.  These  algorithms 
constitute  an  important  subclass  of  Multiple  Access  Algorithms  (MAAs);  they  are 
distributive  and  they  attain  high  throughput  and  low  delays  by  controlling  the  number  of 
simultaneously  transmitting  users.  We  propose  a  method  for  delay  distribution  analysis  of 
Window  RAAs.  It  is  shown  that  a  recently  proposed  method  for  computing  bounds  on  the 
moments  of  the  delay  process-,- -efin  be  extended  to  provide  bounds  on  the  distribution  of  the 
delays.  The  quantities  of  interest  are  related  to  the  solution  of  a  deni”^“rqble  system  of 
lineal,  equations.  The  methodology  is  applied  to  the  delay  distribution  analysis  of  two 
well  known  Window  RAAs. 
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CHAPTER  I 


INTRODUCTION 


1.1  Overview  of  the  Multiple  Access  Problem 

We  consider  the  situation  where  a  number  of  geographically  distributed,  indepen¬ 
dent  users  wish  to  communicate  with  each  other,  or  with  a  central  station,  by  transmitting 
messages  over  a  single  channel. 

The  problem  then  is  the  design  of  an  algorithm,  that  is  a  set  of  rules  which  deter¬ 
mine  the  time  instants  when  a  user  can  transmit  his  messages. 

This  problem  is  referred  to  as  the  Multiple  Access  Problem,  and  the  algorithm 
employed  is  called  a  Multiple  Access  Algorithm  (MAA).  Typical  examples  of  multi  - 
user  communication  systems  are  the  following  [1] : 

(1)  Geosynchronous  Satellite  Systems,  where  many  ground  stations  transmit  to  a 
common  satellite  receiver,  with  the  received  messages  being  relayed  to  the  ground  sta¬ 
tions. 

(2)  Packet  Radio  Networks.  Here  a  message  transmitted  by  a  radio  transmitter 
may  be  received  over  a  wide  area  by  any  number  of  receivers.  This  is  referred  to  as  the 
broadcast  capability.  Thus  a  ground  packet  radio  channel  provides  a  completely  con¬ 
nected  network  topology  for  a  large  number  of  users  within  range  of  each  other. 

It  is  assumed  that  at  most  one  message  can  be  successfully  transmitted  at  a  time 
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simultaneously  over  the  common  channel  we  say  that  a  collision  occurred,  then  the 
received  signal  is  the  sum  of  attenuated  transmitted  signals,  perhaps  corrupted  by  distor¬ 
tion  and  noise.  The  information  contained  in  the  original  messages  is  assumed  lost,  thus 
each  message  involved  in  a  collision  must  be  retransmitted  at  some  later  time,  with 
further  retransmissions  possible  until  the  message  is  successfully  transmitted. 

The  MAA  may  be  centralized  (when  a  central  controller  dictates  the  action  of  each 
user),  distributed  (when  each  user  acts  according  to  prespecified  rules  without  the  inter¬ 
vention  of  a  central  controller),  or  a  combination  of  both. 

The  performance  of  a  MAA  depends  on  the  statistical  nature  of  the  message  gen¬ 
erating  process.  The  Frequency  Division  Multiple  Access  (FDMA) ,  and  Time  Division 
Multiple  Access  (TDMA)  algorithms  are  the  first  designed  to  deal  with  the  problem. 
They  are  both  Collision  -  Free  algorithms. 

In  the  FDMA  algorithm,  the  available  bandwidth  is  divided  in  bands,  and  each  band 
is  dedicated  to  a  user.  In  the  TDMA  algorithm,  the  time  is  divided  in  frames,  each  frame 
is  subdivided  into  slots,  and  each  user  is  allowed  to  transmit  only  once  per  frame,  and 
only  within  the  slot  dedicated  to  him.  When  the  users  Tansmit  messages  continuously, 
both  the  FDMA  and  the  TDMA  are  very  efficient,  permitting  100  %  utilization  of  the 
channel.  However,  if  the  number  of  users  assigned  channels  is  large  and  possibly  time 
varying,  and  if  each  user  transmits  bursts  of  data  infrequently,  both  the  above  algorithms 
induce  unnecessarily  long  delays.  Consider  for  example,  the  TDMA  algorithm  and, 
assume  that  there  are  M  (M»l),  users  in  the  network.  A  message  generated  by  some 
user  will  have  to  wait  an  average  of  M/2  time  slots  for  its  transmission,  even  if  no  other 
user  has  messages  for  transmission. 
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Measurements  conducted  on  time  -  sharing  systems  indicate  that  both  computer  and 
terminal  data  streams  are  bursty,  [2].  Depending  on  channel  speed  the  ratio  between  the 
peak  and  the  average  data  rates  may  be  as  high  as  2000  to  1,  [3].  From  the  above  discus¬ 
sion  it  becomes  clear  that,  if  the  number  of  users  is  large,  and  if  the  users  generate  bursts 
|  of  messages  infrequently,  different  channel  sharing  techniques  are  needed. 

The  Random  Access  Algorithms  (RAAs),  are  designed  exactly  for  this  purpose. 
The  main  difference  between  a  RAA  and  a  perfectly  scheduled  MAA  (  e.g  TDMA,  or 
FDMA)  is  that  the  former  allows  simultaneous  transmission  attempts  by  different  users 
in  the  same  frequency  band.  Therefore,  collisions  are  possible.  Our  objective  then  is  to 
design  Collision  Resolution  Algorithms  (CRAs),  which  resolve  the  resulting  collisions 
and  maintain  stable  system  operation  and  low  delays  (time  that  elapses  from  the  instant 
when  a  message  is  generated  until  this  message  is  successfully  transmitted).  Due  to  the 
fact  that  a  portion  of  the  channel  capacity  is  used  for  the  resolution  of  collisions,  the 
channel  utilization  is  in  general  less  than  100  %.  However,  the  message  delays  at  low 
traffic  rates  are  generally  low,  and  the  channel  utilization  depends  on  the  total  traffic  rate, 
and  not  on  the  number  of  users. 

Only  slotted,  synchronous  Random  Access  Communication  Systems  will  be  con- 
j  sidered.  In  the  sequel  we  list  the  modeling  assumptions  and  briefly  discuss  their  implica¬ 

tions. 

(1)  Slotted,  Synchronous  Systems.  Messages  are  divided  into  constant  size  pack- 
*  ets,  and  each  packet  requires  one  time  slot  for  transmission.  All  users  are  synchronized  to 

the  beginning  of  the  slot  intervals.  Packet  transmission  may  start  only  at  the  beginning  of 
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(2)  Poisson  Arrivals.  It  is  assumed  that  the  cumulative  packet  generating  process  is 
a  Poisson  random  variable,  independent  from  slot  to  slot.  The  Throughput  of  the  RAA  is 
then  defined  as  the  supremum  of  all  Poisson  rates  X,  such  that,  with  probability  one, 
every  packet  is  successfully  transmitted  with  finite  delay. 

(3)  Collision  and  Perfect  Reception.  It  is  assumed  that  if  two  or  more  users 
attempt  to  transmit  their  packets  in  a  given  time  slot,  then  a  collision  occurs  and  all  the 
information  contained  in  the  transmitted  packets  is  destroyed,  thus  retransmission  is  then 
necessary.  If  just  one  user  sends  a  packet  in  a  given  time  slot,  then  the  packet  is  success¬ 
fully  transmitted. 

(4)  Immediate  Feedback  Information.  At  the  end  of  each  slot,  each  user  or 
perhaps  a  subset  of  the  users,  obtain  feedback  information  about  the  transmission  activity 
in  that  slot.  This  feedback  information  is  not  necessarily  error-free. 

(5)  Infinite  Number  of  Users.  The  system  has  an  infinite  number  of  users  and  each 
newly  arriving  packet  arrives  at  a  new  user.  This  assumption  leads  to  a  worst  case  model. 
The  design  of  algorithms  under  this  assumption  is  of  interest,  since  their  performance 
will  not  be  greatly  affected  by  the  number  of  users  in  the  system. 

The  slotted  system  assumption  (1)  leads  to  a  discrete  -  time  system,  thus  simplifying 
the  analysis.  Synchronizing  the  users  is  not  entirely  trivial,  but  can  be  accomplished  with 
relatively  stable  clocks  and  some  guard  time  between  the  end  of  a  packet  transmission 
and  the  beginning  of  the  next  slot.  The  assumption  of  Poisson  arrivals  (2)  approximates 
well  the  case  of  a  large  number  of  independent  bursty  users.  The  assumption  of  Collision 
and  Perfect  Reception  (3)  ignores  the  possibility  of  errors  due  to  noise,  and  also  ignores 
the  possibility  of  Capture  techniques,  by  which  a  receiver  can  capture  one  transmission 


Feedback  (4)  is  quite  unrealistic,  particularly  in  the  case  of  satellite  channels.  For¬ 
tunately  delayed  feedback  although  complicates  random  access  algorithms,  causes  no 
fundamental  problems. 

A  RAA  is  the  combination  of  the  CRA  and  the  First  Time  Transmission  Rule 
(FTTR).  It  is  distributed  in  nature,  since  action  to  be  taken  by  each  user  is  not  based  on 
the  knowledge  of  the  status  of  the  other  users. 

The  oldest  RAA  is  the  ALOHA  [4]  algorithm.  It  was  developed  around  1970  to 
provide  radio  communication  between  the  central  computer  and  various  data  terminals  at 
the  campuses  of  the  University  of  Hawaii.  Each  user  with  a  new  packet  transmits  this 
packet  at  the  beginning  of  the  slot  immediately  after  the  packet  arrival  time.  The  feed¬ 
back  information  is  of  the  acknowledgment  type.  Therefore  only  the  users  that  are 
involved  in  a  collision  are  informed  about  it.  To  resolve  the  collision  the  above  users 
select  random  time  instants  for  the  retransmission  of  their  packets.  As  a  result,  the  system 
contains  two  types  of  users,  at  each  time  slot:  those  that  have  new  packets  to  transmit, 
and  those  whose  packets  have  previously  collided,  and  are  attempting  to  retransmit 
(Blocked  Users).  The  Aloha  algorithm  was  Fust  analyzed  by  Abramson  [4],  The  number 
of  users  was  assumed  infinite,  and  the  number  of  new  packets  generated  by  the  users  at  a 
given  time  slot,  was  assumed  to  be  Poisson  distributed  with  rate  X  packets  per  slot.  How¬ 
ever,  to  simplify  the  analysis,  some  additional  assumptions  were  made  on  the  distribution 
of  the  blocked  users,  and  a  steady  -  state  equilibrium  condition  was  assumed.  Unfor¬ 
tunately,  a  careful  examination  of  the  algorithm  [5]  shows  that  the  system  is  unstable  for 
any  X  >  0,  in  the  sense  that  the  output  rate  eventually  decreases  to  zero. 

The  ALOHA  system  with  a  Fixed  number  M  of  users  has  also  been  studied  [6].  It 


tion  of  the  system  is  ergodic,  if  and  only  if  the  per  packet  retransmission  probability  p 
can  be  chosen  to  satisfy  the  inequality 

X  <  pM(l-p)M_1 

From  this  inequality  the  following  conclusions  can  be  drawn: 

(a)  If  M— then  X  cannot  exceed  e-1 . 

(b)  If  X,  p  and  M  are  such  that  the  above  inequality  is  satisfied,  then  there  exists  Mo, 
such  that  the  system  is  not  stable  for  all  M  £  Mq.  In  other  words  the  ALOHA  algorithm 
is  sensitive  to  the  number  of  users  in  the  system. 

Since  the  ALOHA  algorithm  is  unstable  under  the  infinite  population  model,  the 
natural  question  is  whether  the  design  of  stable  RAAs  is  possible.  The  existence  of  a 
stable  RAA,  for  the  Poisson  packet  generating  process,  was  shown  independently  [7]  and 
[8],  In  those  studies  the  channel  is  assumed  slotted  and  the  feedback  informs  all  the 
users  in  the  system  whether  that  slot  contained  zero,  one,  or  more  than  one  packet 
transmissions  (ternary  feedback).  The  actions  of  the  users  depend  only  on  the  received 
feedback.  It  is  required  that  every  user  observes  the  feedback  channel  independently  if  he 
has  a  packet  to  transmit,  from  the  time  the  system  started  operating  (continuous  feed¬ 
back  sensing).  The  algorithm  can  operate  even  if  the  feedback  is  binary  (informing 
whether  a  slot  contained  more  than  one  packet  transmissions  or  not,  or  in  other  words 
collision  versus  noncollision  feedback),  and  it  can  be  shown  that  it  is  stable  if  X  <  0.345. 
A  modification  of  this  algorithm  increases  its  maximum  throughput  to  0.42.  Under  ter¬ 
nary  feedback,  algorithms  whose  maximum  throughput  is  0.487  [9], [10]  have  been  dev¬ 
ised.  It  can  be  shown  however,  that  in  the  presence  of  feedback  errors  the  latter  algo¬ 
rithms  become  unstable  [11],  while  those  with  the  binary  feedback  are  still  operational  at 
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The  requirement  of  continuous  feedback  sensing  is  usually  unrealistic.  Commonly 
in  practice,  a  user  becomes  active  only  when  he  has  a  packet  to  transmit.  Also  in  the  case 
of  a  failure,  the  continuity  of  channel  sensing  is  interrupted.  It  is  thus  desirable  to  devise 
algorithms  which  require  that  each  user  observes  the  feedback  only  when  he  is  active. 
Those  algorithms  are  called  Limited  Feedback  Sensing  RAAs. 

In  [12]  the  First  such  algorithm  was  presented.  Modifications  that  improve  its  per¬ 
formance  can  be  found  in  [13],  [14]  and  [15].  The  feedback  is  assumed  to  be  either 
binary  or  ternary,  the  algorithms  have  good  delay  characteristics  for  low  arrival  rates  and 
their  maximum  throughput  varies  from  0.36  to  0.4076.  The  analysis  in  [13],  also  shows 
that  the  behavior  of  those  algorithms  in  the  presence  of  feedback  errors  is  quite  satisfac¬ 
tory.  In  [16],  [17]  and  [18]  another  class  of  limited  feedback  sensing  algorithms  was 
presented.  These  algorithms  can  achieve  throughput  as  high  as  the  highest  throughput 
achievable  by  the  known  continuous-feedback-sensing  algorithms,  at  the  expense  of  a 
small  increase  in  delays  and  complication.  They  are  also  not  blocked  (i.e.,  they  do  not 
lead  to  deadlocks),  in  the  presence  of  feedback  errors. 

Three  things  are  worth  noting  here: 

(i)  Both  condnuous  and  limited  feedback  sensing  RAAs  need  more  feedback  sens¬ 
ing  than  that  assumed  by  the  ALOHA  algorithm.  Indeed,  the  ALOHA  algorithm  requires 
that  each  user  know  only  the  outcomes  of  his  own  transmissions.  The  existence  of  stable 
algorithm  for  the  infinite  population  model  subject  to  an  acknowledgement  type  of  feed¬ 
back  is  still  an  open  problem. 

(ii)  The  algorithms  developed  in  [7],  [8],  [9]  and  [10]  are  called  blocked  access 
algorithms.  In  that  case,  all  users  are  able  to  identify  the  end  of  a  Collision  Resolution 
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feedback  continuously,  (a  CRI  is  defined  as  the  time  period  between  the  time  a  collision 
occurred  and  the  time  when  all  the  packets  involved  in  the  collision  have  been  success¬ 
fully  transmitted).  On  the  other  hand  the  algorithms  developed  in  [12],  [13],  [14]  and 
[15]  ,  are  called  free  access  algorithms.  The  latter  algorithms  do  not  need  the  channel 
feedback  history  before  the  packet  generation  to  determine  the  time  of  the  first  transmis¬ 
sion  attempt.  (A  user  with  a  new  packet,  transmits  this  packet  in  the  very  next  slot  fol¬ 
lowing  packet’s  generation  time). 

(iii)  The  algorithms  in  [16],  [17]  and  [18]  exploit  the  fact  that  after  a  user  becomes 
active,  a  finite  number  of  time  slots  is  needed  to  determine  the  status  of  the  channel,  by 
observing  the  feedback. 

For  a  more  detailed  overview  of  the  Random  Access  Problem  and  the  proposed 
algorithms,  the  interested  reader  is  referred  to  [1],  [9],  [11]  and  [19]. 


1.2  OUTLINE  OF  THE  Ph.D  DISSERTATION 


The  work  can  be  divided  in  two  parts.  The  first  part  deals  with  Random  Access  Sys¬ 
tems  where  the  objective  is  not  only  throughput  maximization  and  mean  packet  delay 
minimization. 

In  the  second  part  of  the  work,  a  method  for  delay  distribution  analysis  of  window 
RAAs  is  presented.  In  more  detail,  the  Dissertation  is  organized  as  follows: 

Multiple  Access  Algorithms  are  generally  designed  to  allow  the  successful 
transmission  of  all  the  generated  messages  with  finite  delays.  There  are  applications. 
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threshold. 

In  Chapter  II,  we  propose  and  analyze  a  window  type  RAA  for  time-constrained 
applications.  The  algorithm  guarantees  that  the  total  delay  of  each  successfully 
transmitted  packet  does  not  exceed  a  prespecified  threshold,  and  it  requires  continuous 
feedback  sensing.  The  algorithm  can  be  easily  modified  to  operate  under  limited  feed¬ 
back  sensing.  For  various  thresholds  on  the  transmission  delays,  and  various  input  rates, 
we  computed  the  percentage  of  the  successfully  transmitted  traffic  and  the  expected 
delay  of  a  successfully  transmitted  packet.  We  found  that  significant  improvement  of  the 
expected  packet  delay  is  attained  at  the  expense  of  minimal  traffic  loss. 

In  Chapter  III,  we  consider  the  RAA  proposed  in  Chapter  n  without  imposing  any 
delay  constraints.  We  perform  throughput,  delay,  and  feedback  error  sensitivity  ana¬ 
lyses.  When  the  packet  generating  process  is  Poisson,  the  algorithm  attains  the  same 
throughput  with  that  attained  by  the  Capetanakis’s  Dynamic  Algorithm  [11],  and  it 
induces  lower  mean  packet  delays  for  arrival  rates  above  0.30.  The  algorithm  is  very 
insensitive  to  feedback  errors.  It  operates  in  systems  where  the  packet  generating  process 
is  not  Poisson,  (e.g  where  more  than  one  packets  can  be  generated  at  a  given  time 
instant).  It  is  worth  noting  here,  that  the  algorithms  in  [9],  [10],  are  blocked  when  more 
than  one  packets  can  be  generated  at  the  same  time  instant.  In  the  last  section  of  this 
chapter,  we  analytically  evaluate  the  steady-state  distribution  of  the  distance  between  two 
consecutive,  successful  transmissions,  when  the  proposed  algorithm  is  employed.  The 
latter  distribution  is  important  when  one  studies  interconnected  systems  which  employ 
the  algorithm  for  their  internal  transmissions.  A  typical  example  is  the  multi-hop  prob¬ 
lem,  where  part  of  the  output  (internally  successfully  transmitted  packets)  of  a  random 
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In  Chapter  IV,  a  method  for  delay  distribution  analysis  of  Window  RAAs  is 
presented.  Window  RAAs  constitute  an  important  subclass  of  RAAs;  they  are  distribu¬ 
tive  and  attain  high  throughputs  and  low  delays  by  controlling  the  number  of  simultane¬ 
ously  transmitting  users.  The  throughput  analysis  of  algorithms  in  this  class  is  a  relatively 
easy  task.  The  delay  analysis,  however,  presents  difficulties,  due  mainly  to  the  fact  that 
the  window  size  is  variable  (we  only  impose  a  maximum  window  size),  and  the  compli¬ 
cated  state  space  that  the  algorithms  create.  This  restriction  prohibits  the  application  of 
results  from  standard  queueing  theory  in  the  delay  analysis. 

We  show  that  the  methodology  employed  in  [20],  can  be  extended  to  provide 
bounds  on  the  distribution  of  the  delays.  The  quantities  of  interest  are  related  to  the 
solution  of  a  denumerable  system  of  linear  equations.  Methods  for  the  computation  of 
the  constant  terms  and  the  coefficients  of  the  unknowns  of  the  system  are  developed.  The 
methodology  is  applied  to  the  delay  distribution  analysis  of  the  Capetanakis’s  Dynamic 
algorithm  [11],  and  the  Part-and-Try  algorithm  [9]  and  [10],  both  under  binary  (collision 
versus  noncollision)  feedback.  It  can  also  be  applied  directly  to  other  Window  RAAs 
with  different  feedback.  An  interesting  result  of  the  analysis  is  that  as  the  arrival  rate 
increases,  the  tails  of  the  delay  distribution  become  longer,  but  the  median  grows  much 
slower  than  the  expected  delay. 

In  the  last  Chapter,  Chapter  V,  we  consider  a  two-cluster  packet  radio  network.  In 
each  cluster  a  single  forward  channel  is  available  for  packet  transmission.  The  frequency 
bands  dedicated  to  each  cluster  are  different.  Consequently,  a  packet  transmission  over 
the  forward  channel  in  one  cluster  does  not  interfere  with  a  packet  transmission  over  the 
forward  channel  of  the  other  cluster.  Each  cluster  contains  local  and  marginal  users.  The 
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The  marginal  users  are  located  in  the  overlapping  region  of  the  two  clusters.  They  are 
capable  of  transmitting  their  packets  over  either  one  of  the  forward  channels  and,  receiv¬ 
ing  feedback  information  from  either  one  of  the  feedback  channels.  Due  to  the  double 
exposure  the  marginal  users  have  a  choice:  they  can  join  either  cluster  1  or  cluster  2  for 
the  transmission  of  their  packets.  In  this  Chapter,  we  propose  a  dynamic  protocol  accord¬ 
ing  to  which  a  marginal  user  can  select  the  forward  channel  over  which  he  will  transmit 
his  packet.  The  protocol  requires  that  each  marginal  user  with  a  packet  for  transmission 
observes  both  the  feedback  channels  from  the  packet’s  generation  time.  No  a  priori 
knowledge  of  the  input  traffic  rates  or  the  state  of  the  system  is  required. 
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CHAPTER  H 


A  CONTINUOUS  SENSING  WINDOW  RANDOM-ACCESS  ALGORITHM  FOR  MESSAGES 

WITH  STRICT  DELAY  CONSTRAINTS 


D.1  INTRODUCTION 

We  consider  the  case  where  a  number  of  geographically  distributed  stations  wish  to 
communicate  by  transmitting  messages  over  a  single  channel.  The  problem  then  is  the 
design  of  efficient  algorithms  for  allocating  the  channel  among  the  users.  This  problem  is 
referred  to  as  the  Multiple-Access  problem,  and  the  algorithm  employed  is  called  a 
Multiple- Access  Algorithm  (MAA).  The  main  performance  measures  of  an  algorithm  are 
generally  throughput,  induced  delays,  and  error  sensitivity.  In  addition,  Multiple-Access 
algorithms  are  generally  designed  to  allow  the  successful  transmission  of  all  the  mes¬ 
sages  generated  in  the  system  with  finite  delays. 

There  are  applications,  however,  where  a  message  is  considered  lost  if  its  delay 
exceeds  a  certain  threshold.  Typical  examples  of  such  applications  are  the  following: 

a)  The  transmission  of  control  messages  that  are  required  to  monitor  remote  devices 
within  a  specified  amount  of  time. 

b)  Distributed  target-tracking,  where  a  user  attempts  to  track  a  moving  target  by  using 
local  observations  and  the  information  received  from  the  other  users. 

_  \  .  •  •  J*  •  *  •  _ «•_  _ • _ 9 _ _ ■ - . - 
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Recently,  there  has  been  considerable  interest  in  the  design  and  analysis  of  algo¬ 
rithms  suitable  for  time-constrained  applications.  Those  algorithms  disregard  the 
transmission  of  a  message  if  its  delay  exceeds  the  given  threshold.  Therefore,  a  percen¬ 
tage  of  the  generated  traffic  is  lost. 

Random  Access  Algorithms  constitute  an  important  subclass  of  Multiple  Access 
Algorithms,  and  it  is  of  considerable  interest  to  design  such  algorithms  that  will  operate 
efficiently  and  reliably  under  strict  constraints  on  the  delays  of  the  generated  messages. 
As  discussed  in  [21],  the  advantages  of  the  Random  Access  Algorithms  for  time  con¬ 
strained  communications  are  "guaranteed  connectivity,  the  transfer  of  temporary  over¬ 
loading  into  packet  delay  rather  than  blocked  connections,  and  the  ability  to  trade  mes¬ 
sage  loss  for  message  delay".  The  major  disadvantage  is  that  a  Random  Access  Algo¬ 
rithm  induces  variable  message  delay. 

Simulation  and  experimental  studies  of  the  performance  of  the  Ethernet  Protocol 
under  time  constraints  can  be  found  in  [32],  [33],  Window  Access  protocols  for  time  con¬ 
strained  communication  have  been  considered  in  [22],  In  [22],  it  was  assumed  that  con¬ 
straints  are  imposed  only  on  the  time  elapsed  from  the  generation  time  of  a  message  to 
the  time  of  its  first  transmission  attempt.  No  restriction  was  imposed  on  the  additional 
delays  experienced  by  a  message  when  its  first  transmission  attempt  is  not  successful. 
Furthermore,  an  approximate  model  of  the  system  was  developed  and  analyzed.  This 
model,  however,  fails  to  take  into  account  the  fact  that  window  sizes  are  necessarily  vari¬ 
able  in  a  real  system. 

In  this  chapter,  we  propose  and  analyze  a  window  Random-Access  algorithm  for 
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that  the  users  observe  the  channel  feedback  continuously  (continuous  feedback  sensing). 
We  point  out  that  the  algorithm  can  be  easily  modified  to  operate  under  limited  feedback 
sensing,  where  the  users  are  required  to  observe  the  feedback  only  whenever  they  have  a 
message  to  transmit.  However,  the  analysis  of  the  algorithm  is  much  more  involved  in 
the  latter  case. 

This  chapter  is  organized  as  follows:  In  section  n.2,  we  present  the  system  model. 
In  .section  II.3,  we  present  the  operation  of  the  algorithm.  In  section  n.4,  we  extend  the 
basic  ideas  presented  in  [20]  to  provide  an  exact  analysis  of  the  time-constrained  algo¬ 
rithm.  The  performance  of  the  algorithm  is  discussed  in  section  II. 5. 


11.2  SYSTEM  MODEL 

The  transmission  delay  of  a  message  in  a  communication  network  is  defined  as  the 
time  period  between  the  instant  that  the  message  is  generated  and  the  instant  when  its 
successful  transmission  has  been  completed.  When  a  strict  upper  bound  on  the  transmis¬ 
sion  delays  is  imposed,  some  messages  are  lost.  Important  performance  criteria  in  this 
case  are  the  fraction  of  the  lost  message  traffic  and  the  expected  transmission  delays  of 
the  transmitted  messages. 

In  this  chapter, we  consider  a  random  access  packet  network  with  upper  bounds  on 
the  packet  delays.  In  particular,  we  assume  that  a  common  slotted  channel  is  shared  by 
packet  transmitting  independent  users,  that  the  feedback  per  slot  is  binary  CNC  (collision 
versus  noncollision),  and  that  the  feedback  channel  is  errorless.  We  also  assume  that  col¬ 
lisions  result  in  total  loss  of  all  the  involved  packets,  a  packet  transmission  may  start  only 
at  the  beginning  of  some  slot,  and  that  there  are  no  propagation  delays.  Finally,  it  is 
assumed  that  the  packet  generating  process  is  Poisson.  This  is  a  reasonable  assumption 
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since  we  assume  a  large  number  of  independent  bursty  users.  Furthermore  in  [28]  we 
prove  that  for  a  large  class  of  Random  Access  algorithms  (RAAs),  and  given  a  fixed 
value  for  the  cumulative  packet  generation  rate,  as  the  user  population  increases  the  sta¬ 
bility  of  an  algorithm  in  the  class  is  determined  by  its  throughput  under  the  Poisson  user 
model.  For  the  above  system  model  we  propose  and  analyze  a  continuous  feedback  sens¬ 
ing  random  access  algorithm.  Additionally  the  method  we  develop  is  also  applicable  to 
the  case  where  the  number  of  packets  generated  per  slot  forms  an  i.i.d  process  with  gen¬ 
eral  distribution  (this  may  be  the  case  in  some  synchronized  systems),  and  to  the  finite- 
population  user  model. 


II.3  THE  ALGORITHM 

Let  T  be  the  upper  bound  on  the  packet  transmission  delay,  and  let  it  be  common  to 
all  users.  Let  t  be  a  time  instant  corresponding  to  the  beginning  of  some  slot,  such  that, 
for  some  tx  <  r,  all  the  packet  arrivals  in  (0,  rj  ]  have  been  either  successfully  transmitted 
or  rejected  by  the  algorithm,  and  there  is  no  information  regarding  the  arrival  interval 
(r | ,  r]  ,  (Fig.II.l).  Such  a  point  t  is  called  a  "collision  resolution  point”  (CRP).  Let  time 
be  measured  in  slot  units.  Then,  in  slot  t,  the  arrivals  in  (t 2,  13]  transmit,  where  for  A 
being  an  algorithmic  constant  we  have:  1 2  =  max  (t  -  (T- 1),  t ] )  and  1 3  =  min(r2  +  A,  r). 
At  the  same  time,  the  arrivals  in  (t)  ,12]  are  automatically  rejected,  because  their 
transmission  delay  has  exceeded  T.  The  arrival  interval  (0,  r  1  ]  above  is  called  "resolved 
interval,"  the  interval  (f2,  I]  is  called  the  "lag  at  t",  and  the  arrival  interval  (r2-  r 3]  is 
called  the  "examined  interval."  If  (f2,  13]  contains  zero  or  one  packets,  then  it  is 
resolved  at  slot  t.  If  (f2,  *3]  contains,  instead,  at  least  two  packets,  then  a  collision 
occurs  at  t,  whose  resolution  begins  at  slot  t  +  1.  Until  the  collision  at  t  is  resolved,  no 
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the  resolution  of  the  latter  collision  is  called  the  "collision  resolution  interval"  (CRI).  Let 
d=t-t2  be  the  length  of  the  lag  at  t.  The  length  of  the  collision  resolution  interval  is  not 
allowed  to  exceed  A  =T -d.  Thus,  packets  in  O2,  f 3],  that  are  not  successfully  transmitted 
by  time  t+A,  are  rejected.  This  last  rule  assures  that  all  the  successfully  transmitted  pack¬ 
ets  have  delay  less  than  T.  Note  that  under  this  rule  it  is  possible  that  a  packet  (e.g., 
packet  2  in  Fig.  II.  1)  with  delay  T- A  ,  and  not  T,  be  rejected.  However,  as  we  shall  see  in 
section  5,  the  optimum  value  of  A  is  close  to  2  slots  and,  therefore,  this  rule  does  not 
represent  any  severe  restriction  in  practice. 

Let  x ,  denote  the  feedback  corresponding  to  slot  t,  where  either  xt  =  c,  for  collision 
slot,  or  x,  =  nc  ,  for  noncollision  slot.  During  some  CRI,  each  involved  user  acts 
independently,  utilizing  a  counter,  whose  value  at  time  t  is  denoted  by  rt.  When  a  user 
transmits  a  packet  for  the  first  time  he  sets  r,  =  1 .  The  value  of  r,  is  updated  and  used  as 
follows: 

1.  The  user  transmits  each  time  t  that  r,  =  1  .  The  corresponding  packet  is  successfully 
transmitted  at  t  if  and  only  if, 

r,  =  1  and  xt  =  nc 

2.  The  transitions  in  time  of  the  counter  value  r,  are  as  follows: 

(a)  If  .r,_i  =  nc  and  =  2  ,  then  r,  =  1. 

(b)  If  x,-\  =  c  and  r,_!  =  2,  then  r,  =  2. 

(c)  If  xt-\  =  c  and  r,_\  =  1,  then 


1 


1  ;  w.p.  0.5 

2  ;  w.p.  0.5 
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The  algorithmic  operation  can  be  depicted  by  a  two  cell  stack.  At  each  slot  r ,  cell  1 
contains  the  transmitting  users  (those  with  r,  =  1  ),  and  cell  2  contains  the  withholding 
users  (those  with  r,  =  2  ).  In  contrast  Capetanakis’s  algorithm  distributes  the  unsuccessful 
users  across  the  cells  of  an  infinite  cell  stack.  The  algorithm  allows  every  user  in  the  sys¬ 
tem  to  determine  accurately  the  time  instant  t'  when  a  CRI  has  ended;  a  CRI  ends  either 
when  two  consecutive  noncollision  slots  occur,  or  when  its  length  equals  T-d  where  d  is 
the  length  of  the  lag  at  the  beginning  of  the  CRI.  Then,  the  algorithm  reinitializes  with  a 
new  examined  interval.  We  note  that  the  algorithmic  parameter  A  is  subject  to  optimal 
selection,  for  performance  enhancement. 


1 

HM - -  d  <  T  - 


1 


lu,d<T-d 


Figure  II.  1 

Illustration  of  the  relationships  among  certain  random 
variables  related  to  the  operation  of  the  algorithm. 
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II.4  ANALYSIS 

Consider  the  algorithm  in  section  n.3.  Let  the  system  start  operating  at  time  zero, 
and  let  us  consider  the  sequence  in  time  of  lags  that  are  induced  by  the  algorithm.  Let  Xt 
denote  the  length  of  the  i-th  lag,  where  i>l.  Then,  the  first  lag  corresponds  to  the  empty 
slot  zero;  thus,  Xj  =  1.  In  addition,  the  sequence  X,;  iSl  is  a  Markov  chain  whose  state 
space  is  at  most  countable.  Let  Dn  denote  the  delay  experienced  by  the  n-th  successfully 
transmitted  packet  arrival,  as  induced  by  the  algorithm;  that  is,  the  time  between  the 
arrival  of  the  packet  and  its  successful  transmission.  Let  the  sequence  Tit  it  l  be  defined 
as  follows:  Let  T\  =  1  ,  and  define  7i  + j  as  the  first  CRP  after  7,  at  which  the  lag  has 
length  one. 

Let  Rj,  it  1  and  F\,  it  1  denote  respectively  the  number  of  successfully  transmitted 
packets  and  the  number  of  rejected  packets  in  the  time  interval  (0, 7,].  Then, 
Qi  =Ri+ 1  ~Ri *2:1  and  G,  -7,,  j>1  denote  respectively  the  number  of  success¬ 

fully  transmitted  and  the  number  of  rejected  packets  in  the  interval  (7,  ,  7l+i],  The 
sequences  Q,,  it  1  and  G,,  <>1  are  sequences  of  i.i.d.  random  variables;  thus  Rt,  it l  and 
7j,  />1  are  renewal  processes.  In  addition,  the  delay  process  £>„,  n>l  induced  by  the 
algorithm  is  regenerative  with  respect  to  the  process  7,,  it  1  and  the  distribution  of  (2,  is 
nonperiodic,  since  P  (Qi  =  1)  >  0. 

Let  us  define, 

Q^ 

Z  =  E{Q\}  ,  W  =  E{^Dt),  H  =  E{T2-Tx }  (II.l) 

j=i 

From  the  regenerative  arguments  in  [201,  it  follows  that  the  fraction,  p,  of  success¬ 
fully  transmitted  packets,  and  the  expected  steady-state  delay,  D,  of  the  successfully 
transmitted  packet  are  respectively  given  by  the  following  expressions: 


r 


- 
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p  =  za//r1 

(II.2) 

t 

D  =  WZ~l 

(II.3) 

where  X  is  the  intensity  of  the  Poisson  packet  generating  process. 

Towards  the  computation  of  the  expected  values  Z,  H,  and  W,  let  us  define  the  fol¬ 
lowing  quantities  (see  Figure  II.  1),  where  t2,  t2,  and  t  are  as  in  section  H3: 

nuj  :  Number  of  packet  arrivals  in  [t2  ,  f 3 )  that  are  successfully  transmitted 

during  the  collision  resolution  process,  given  that  t2-t2  =  m.,  and  t-t2=d 

Sum  of  the  delays  of  the  nuj  packets,  after  time  t. 

Sum  of  the  delays  of  the  nu  d  packets,  until  the  instant  r3  . 

The  number  of  slots  needed  to  examine  an  interval  of  length  u  given  that 
t-t2=d.  Note  that  lUtd  <  T-d. 

Conditional  expectation  of  the  random  variable  X,  given  that  the  length  of 
the  initially  transmitted  interval  is  u,  and  t-t2=d. 

hd  :  The  number  of  slots  needed  to  reach  a  CRP  with  lag  equal  to  one.when 

starting  from  a  CRP  with  lag  equal  to  d  . 

wd  :  The  cumulative  delay  experienced  by  all  the  packets  that  were  success¬ 

fully  transmitted  during  the  hd  slots. 

ad  :  The  number  of  packets  that  are  successfully  transmitted  within  the  interval 


?u,d  ■ 

Vu.d  : 

^u.d  ■ 

E{X  \u,d}  : 
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P(l  I  u,d ) :  Given  that  the  interval  to  be  examined  has  length  u  and  t  - t2=d l,  the  proba¬ 
bility  that  the  corresponding  collision  resolution  interval  has  length  /. 

Hd  =  E{hd} 

Wd  =  E{wd}  (II. 4) 


Ad  =  E{  ad} 

We  note  that  the  quantities  in  (2),  (3),  and  (4)  are  such  that, 
Z  =  A ! ,  H  =  H  i ,  and  W  =  W  j .  Denoting  by  xd  either  one  of  the  random  variables 
hd ,  wd,  ad,  the  operations  of  the  algorithm  in  section  II. 3  induce  the  following  relation¬ 
ships. 

-\T-d\>ldd>\ 

Xd  =  +  xd.M^\ T~\>d>A 

;  where  [  a\  denotes  integer  part  of  a ,  and 


l minCA.d).d  1  for  the  r.v.  hd 

®d  =1  ¥mm <b.d),d  +  +  m^d-A^^n^du  ;  for  the  r.v.  wd  (II. 6) 

«miiy'A,d;.d  ;  for  the  r.v. ad 

Taking  expectations  in  (5),  and  denoting  X^Efx^,  we  obtain: 

L  t-4 

Xd  =  E{Qd}+  ZXmP(m\d,d)  ;  1  <d<A 

m=2 


(H.7) 

[T-4 

Xd  =  E{ed}  +  2  Xd-6+m  P(m  I A ,d)  :  A  <  d  <  T-l 

m=\ 

;  where  d  takes  at  most  denumerable  values  in  [1,7'-1],  In  (7),  the  index  m  takes 


values  up  to  [T-d]  because  lUtd  S  T-d. 


21 


I 


I 
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Note  that  the  system  in  (7)  is  a  finite  system.  Moreover,  if  A  is  an  integer,  then  d 
can  only  take  integer  values  in  [1,  T- 1].  In  either  case,  the  system  is  of  the  following 
form: 

x  =  b  +  @x  or  (IS) x  =  b 

where  0  is  a  strictly  substochastic  irreducible  matrix.  It  follows  [34]  that  if  b> 0,  then  the 
system  in  (7)  has  a  unique  nonnegative  solution.  In  Appendix  A,  we  include  the  recur¬ 
sions  pertinent  to  the  computation  of  E{Qdi  and  P(l  \u,d).  The  latter  recursions,  in  con¬ 
junction  with  the  system  in  (7),  provide  the  fraction  of  the  successfully  transmitted  traffic 
and  the  expected  delay  of  the  successfully  transmitted  packets,  for  given  T,  A  and  X  . 
Given  the  latter  parameters,  the  quantities  in  (2)  and  (3)  are  functions  of  them;  we  will 
denote  them  by,  pf(A,X)  and  Dj-fAAj,  respectively. 


II.5  PERFORMANCE  EVALUATION 


The  value  of  A  is  a  basic  design  parameter  for  window  Random  Access  Algorithms. 
If  no  delay  constraints  are  imposed,  A  is  chosen  so  that  the  stability  region  of  the  algo¬ 
rithm  is  maximized.  For  systems  with  strict  delay  constraints,  an  upper  bound  T  on  the 
transmission  delay  of  each  successfully  transmitted  packet  is  given.  Since  a  fraction  of 
the  generated  packets  is  then  rejected,  an  important  quality  measure  in  this  case  is  a 
lower  bound,  ej ,  on  the  fraction  of  the  successfully  transmitted  traffic.  Given  T and  e  \  as 
above,  the  value  of  A  is  determined  from  the  following  optimization  problem  :  Deter¬ 
mine  A,  so  that  the  input  traffic  rate  is  maximized  while  the  fraction  of  the  successfully 
transmitted  traffic  remains  greater  than  e  t .  That  is,  the  quantity  Xf  below  is  sought. 


Xt  .  = 


'  Or  <  A  1  V  >  0  .  t 


Kurt 
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Due  to  the  complexity  of  the  expression  that  provides  the  system  parameter 
Pr(A,  X),  it  is  hard  to  solve  the  optimization  problem  in  (8)  in  an  analytical  fashion.  For 
fixed  A,  however,  it  is  simple  to  compute  the  quantity  Xf,  (A)  defined  below. 

(A)  =  sup  (X  :  Pt-(A.X)  >  ej)  (II.9) 

In  table  n.  1  we  exhibit  the  X^  ei(A)  values,  for  A=  1.5,  2,  2.5,  3,  and  various 

values  of  the  parameters  T  and  e  i .  From  the  above  table,  we  observe  that  the  algorithm  is 
relatively  insensitive  to  the  selection  of  the  parameter  A,  for  the  A  e  [2,  3]. 

If  in  addition  to  e  i ,  it  is  also  desirable  to  impose  an  upper  bound  e2  on  the  expected 
delay  of  the  successfully  transmitted  packet,  then  the  value  of  A  is  determined  from  the 
following  optimization  problem:  Select  A  so  that  the  quantity  Xft  ev  defined  below,  is 
attained. 

^7\  <»i.  (X  :  Pr(A,  X)  >  e\  ,  D^fA,  X)  <  e2)  (11.10) 

The  analytic  determination  of  X^  is  again  not  feasible.  We  thus  compute  instead 
the  rates  Xf  ej(A)  defined  below. 

<r,.<2(A)=  sup  (\:pT(A,\)>e],DT(A,\)<e2)  (II.  11) 

In  table  H.2  we  exhibit  Xf  ej(A)  values,  for  A  =  1.5,  2,  2.5,  and  3,  and  for  various 

values  of  the  parameters  T,  ej  and  e2-  For  A  €  [2,3],  we  observe  again  insensitivity  of 
the  algorithm  to  the  selection  of  the  parameter  A. 

From  table  II.  1  we  observe  that  for  fixed  e  j  and  A,  Xf  ei  (A)  is  an  increasing  func¬ 
tion  of  T.  In  contrast,  from  table  II. 2  we  observe  that  for  fixed  eL  e2t  and  A,  there  exists 
some  Tmaxfe!,  e2  A),  such  that  Xf  ei  f2(A)  increases  for  1  <T  £Tmax(e i,  e2.  A)  and 
decreases  for  T>  TmAX(e\t  e2  A).  The  latter  observation  can  be  explained  as  follows: 
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cessfully  transmitted  traffic.  For  large  T  values,  on  the  other  hand,  the  determining  fac¬ 
tor  is  the  constraint  on  the  expected  delay  of  the  successfully  transmitted  packet. 

It  is  interesting  to  compare  the  performance  of  the  algorithm  under  strict  delay  con¬ 
straints,  to  its  performance  when  no  such  constraints  are  imposed.  In  the  latter  case,  for 
input  traffic  rates  within  the  stability  region  of  the  algorithm,  the  fraction  p  of  the  suc¬ 
cessfully  transmitted  traffic  equals  one,  and  the  throughput  represents  then  the  maximum 
maintainable  rate  of  the  input  traffic.  In  Appendix  A,  we  compute  relevant  useful 
bounds.  Using  those  bounds  and  the  methodology  in  [20],  we  compute  the  throughput 
p*,  the  optimal  parameter  A*  that  attains  p*,  and  the  expected  delays  D„(A*,  A.)  for 
te(0,p*).  We  found: 

p*  =p*(A*)  =  0.429  ,  A*  =2.33  (11.12) 

Figure  n.2,  illustrates  expected  delays  for  various  values  of  T  and  for  A=2..  In  the 

same  Figure  we  illustrate  the  delays  induced  when  no  delay  constraints  are  imposed.  In 
Figures  II. 3  and  n.4  we  illustrate  the  fraction  of  the  successfully  transmitted  traffic  and 
the  successfully  transmitted  traffic  respectively,  for  various  values  of  T  and  for  A=2.  We 
observe  that  for  input  traffic  rates  above  0.30,  imposing  strict  constraints  on  the  delays  of 
the  successfully  transmitted  traffic  results  in  significant  improvement  on  the  expected 
delay  of  the  successfully  transmitted  packet;  the  price  paid  then  is  loss  of  some  of  the 
generated  traffic.  For  example,  for  X.  =  0.4,  A  =  2,  and  7=20,  the  algorithm  attains 
P20U,  0.4)  »  0.95  ,  and  D20  (2,0.4)  =  7.6,  versus  D»  (A*'  0.4)  =>  19  in  the  absence  of 
an  upper  bound  T,  while  for  T  =  20  the  delay  of  each  successfully  transmitted  packet  is 
simultaneously  not  exceeding  20  slots.  The  tradeoff  between  acceptable  delays  and 
acceptable  loss  in  traffic  is  determined  by  each  particular  application. 
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feedback  [11].  Compared  to  Capatanakis’s  algorithm,  the  proposed  algorithm  has  the 
property  that  both  the  unconstrained  and  the  constrained  version  can  be  easily  modified 
to  operate  under  limited  feedback  sensing.  This  is  because  during  a  collision  resolution 
process,  no  more  than  two  consecutive  noncollision  slots  can  be  encountered.  This  pro¬ 
perty  can  be  used  to  provide  a  limited  feedback  sensing  version  of  the  algorithm 
presented  in  this  chapter  (see  also  [17],  [18]).  It  is  also  important  to  note  that  due  to  the 
very  simple  rules  of  the  collision  resolution  process,  the  algorithm  is  less  sensitive  to 
feedback  channel  errors  than  the  Capetanakis’s  dynamic  algorithm.  A  complete  analysis 
of  the  unconstrained  version  of  the  algorithm  is  presented  in  Chapter  III.  The  following 
argument  can  be  used,  however,  to  justify  our  claim:  Let  5  be  the  probability  that  an 
empty  slot  is  erroneously  interpreted  as  a  collision  slot,  and  assume  that  this  is  the  only 
form  of  error  on  the  feedback  channel.  Let  B^\  E^\  be  respectively  the  average 
number  of  collision,  empty  and  successful  slots  during  a  multiplicity  k  collision  resolu¬ 
tion  process.  Let  L^5)  be  the  average  number  of  slots  needed  for  a  multiplicity  k  collision 
to  be  resolved.  It  is  easy  to  see  that 

fii8)=fii0),  Si6)=Si°\  E\ fc8)  =  40)7r^r  ;  0  <  5  <  1,  k>2  (11.13) 

(l-o) 

It  follows  that 

Li8)  =  fif  +S<k0)+Ei0)  y^y  =  L^+Ep  yy^y  <  yj^y  ;  0<5<  1 ,  k>2  (II.  14) 
It  can  also  be  seen  that 

45)=— tt  :  0  <  5  <  1 
(1-5)2 

oo 

If  fhM-  X  L<f)e*xklk !,  the  throughput  of  the  algorithm  is 

*= o 


I 


(11.15) 
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Xs>Xo(l-5)2  (11.16) 

Therefore,  the  throughput  of  the  algorithm  remains  positive  for  0  <  5  <  1.  It  is  well 

known  [11],  however,  that  the  throughput  of  Capetanakis’s  algorithm  reduces  to  zero  for 

5  >0.5. 

Gallager’s  algorithm  [11],  can  also  be  modified  to  operate  under  delay  constraints. 
The  method  developed  in  this  chapter  can  be  applied  to  the  analysis  of  the  latter  algo¬ 
rithm,  but  the  computations  become  much  more  involved  and  greater  care  is  needed 
since  the  developed  systems  of  linear  equations  become  infinite.  For  systems  where  the 
Poisson  user  model  is  valid,  it  is  expected  that  the  modification  of  Gallager’s  algorithm 
will  provide  improved  performance.  One  advantage  of  the  algorithm  presented  here,  is 
that  it  can  operate  in  systems  where  the  Poisson  user  model  is  not  valid  (e.g.  when  more 
than  one  packets  are  generated  at  a  specific  instance  due  to  synchronous  transmission),  in 
contrast  to  the  algorithm  in  [11],  and  that  analytical  evaluation  of  its  performance  and 
optimal  choice  of  the  design  parameter  A  is  then  feasible. 

n.6  CONCLUSIONS 

We  considered  the  case  where  strict  constraints  on  the  delays  of  successfully 
transmitted  packets  are  imposed.  For  this  case,  we  presented,  analyzed,  and  evaluated  a 
simple  random  access  algorithm.  The  algorithm  requires  continuous  feedback  sensing. 
We  analyzed  the  algorithm  under  the  Poisson  user  model  assumption.  The  analysis  of  the 
algorithm  is  based  on  its  regenerative  characteristics.  For  various  bounds  on  the 
transmission  delays  and  for  various  input  rates,  we  computed  the  percentage  of  the  suc¬ 
cessfully  transmitted  traffic  and  the  expected  delay  of  the  successfully  transmitted 
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CHAPTER  III 


A  WINDOW  RANDOM  ACCESS  ALGORITHM 
FOR  MULTI-USER  PACKET  RADIO  SYSTEMS 


IH.1  INTRODUCTION 

In  data  networks  with  bursty  users,  the  most  appropriate  multiple-access  techniques 
are  those  belonging  to  the  class  of  Random  Access  Algorithms  (RAAs).  This  is  due  to 
the  fact  that  RAAs  are  distributed,  and  they  induce  low  delays  for  low  input  rates. 

The  important  performance  characteristics  of  a  RAA  are:  throughput,  induced 
delays,  and  sensitivity  to  feedback  channel  errors.  Insensitivity  to  feedback  errors  is 
important  when  operating  under  noisy  conditions.  Different  sources  of  noise  include:  (i) 
interference  between  closely  located  cells,  which  use  the  same  frequency  band  (multi  hop 
networks),  or,  (ii)  multi-path  fading.  The  effects  of  noise  can  be  modeled  as  erroneously 
observed  feedback,  see  [35]. 

In  this  chapter,  we  consider  the  RAA  presented  in  chapter  II  without  imposing  any 
delay  constraints;  i.e.,  T  =  We  then  perform  throughput,  delay,  and  feedback  error 
sensitivity  analyses.  We  also  evaluate  the  output  traffic  interdeparture  distribution 
induced  by  the  algorithm. 

Under  the  Poisson  user  model  assumption,  the  algorithm  attains  the  same 
throughput  with  that  attained  by  the  Capetanakis’s  dynamic  algorithm,  [11],  while  it 
induces  lower  delays  for  arrival  rates  above  0.30.  It  also  exhibits  better  resistance  to  feed- 
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back  errors.  The  evaluation  of  the  output  traffic  interdeparture  distribution  is  important, 
in  the  case  of  interconnected  networks  each  of  which  employs  the  RAA  for  internal 
transmissions.  The  evaluation  of  the  latter  distribution  when  either  the  Capetanakis’s, 
[1 1],  or,  the  Gallager’s,  [11],  algorithms  are  used  remains  an  open  problem. 

The  organization  of  this  chapter  is  as  follows:  In  section  III.2,  we  present 
throughput  and  delay  analysis.  In  section  III.3,  we  present  feedback  error  sensitivity 
analysis,  and  discuss  the  operation  of  the  algorithm  under  limited  feedback  sensing.  In 
section  III.4,  the  output  traffic  interdeparture  distribution,  induced  by  the  algorithm,  is 
evaluated.  In  section  III.5,  some  conclusions  are  drawn. 

in.2  THROUGHPUT  AND  DELAY  ANALYSES 

In  this  section,  we  present  throughput  and  delay  analyses  of  the  algorithm.  We 
assume  continuous  feedback  sensing  and  the  Poisson  user  model. 

Consider  the  algorithm  in  Section  II.3.  Since  we  do  not  impose  any  delay  con¬ 
straints,  we  have  that  t\  =ti .  Let  the  system  start  operating  at  time  zero.  Let  r,  ;  i  >  1  be 
the  sequence  of  successive  CRPs,  and  Xj  be  the  lag  at  f;.  The  sequence  X,  ;  i  >  1  is  a 
Markov  Chain  with  state  space  F,  F  is  an  at  most  denumerable  subset  of  the  interval  [1, 
«> ).  It  can  be  seen  that  any  state  can  be  reached  from  any  other;  therefore  X;  ;  i  >  1  is  an 
irreducible  Markov  Chain.  Since  P  (Xi+i  =  1 1  X,  =  1)  >  0  ,  we  conclude  that  X,;  i  >  1  is 
also  aperiodic.  Therefore  Pake’s  Lemma  [30]  applies,  and  gives  that  the  following  condi¬ 
tion  is  sufficient  for  the  ergodicity  of  the  Markov  Chain,  (stability  of  the  system): 

E  (/  I  A,cf)  <  A  (ffl.1) 

where  E  (/  I  A ,d)  denotes  the  expected  length  of  a  CRI,  given  that  it  starts  with  an  exam- 

nf  tonath  A  and  with  a  lag  of  length  d. 
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Since  the  Markov  Chain  is  uniformly  downward  bounded  (there  exists  a  constant  m 
such  that  the  transition  probabilities  p kj  satisfy  p kj  =  0  for  j  <  k  -  m.  Here  m  =  A-  1), 
Kaplan’s  Theorem  [31]  applies  and  gives  that: 

If 

£(/  \A,d)  >  A 

then  the  Markov  Chain  is  not  ergodic  and  the  system  is  unstable. 

Let  Lk  denote  the  expected  length  of  a  CRI  given  that  it  starts  with  a  collision  of 
multiplicity  k.  We  can  then  write: 

E{1  \A,d)  =  ££(/  IA ,d,k)e~^^-=  (HI.2) 

k= 0  K-  *  = 0  K- 

since 

£(/  \A,d,k)  =  Lk  (III. 2a) 

depends  only  on  k. 

In  Appendix  B  we  show  that: 

(i)  Lk  ■,  k  >0  can  be  computed  recursively,  and 

(ii)  Lk  are  quadratically  upper  bounded,  Lk  <  Lk  =  cotz  +  $k  +  y ;  k  >  2. 

Expression  (1)  together  with  (i)  and  (ii)  are  used  in  the  computation  of  the  algorithmic 
throughput.  The  details  of  the  analysis  are  presented  in  Appendix  B. 

We  define  the  delay  Dn,  experienced  by  the  n-th  packet  as  the  time  difference 
between  its  arrival  and  the  end  of  its  successful  transmission.  We  are  interested  in 
evaluating  the  first  moment  of  the  steady  state  delay  process,  when  it  exists.  Let  Tt  =  1  , 
X !  =  1  ,  and  define  T,+i  as  the  first  CRP  after  7,  at  which  the  lag  has  length  one.  From 
the  description  of  the  algorithm  it  can  be  seen  that  the  induced  delay  process  probabilisti¬ 
cally  restarts  itself  at  the  beginning  of  each  slot  T,  ,  i  =  1  ,2  ,  ...  The  interval  [7,,Tl+i) 
will  be  referred  to  as  the  i-th  session.  Note  that  the  sessions  have  lengths  that  are  i.i.d  ran- 
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dom  variables. 

Let  Ri  ;  i  =  1 ,2 ,  ..  ,  denote  the  number  of  packets  successfully  transmitted  in  tne 
interval  (0,  Tt  ]  ;  (note  that  R,  also  represents  the  number  of  packets  arrived  during  the 
interval  [0,  7",  -  1),  since  7",  is  a  CRP  at  which  the  lag  has  length  one).  Then, 
Qi  =/?l+1  -  Ri ;  i  >  1,  is  the  number  of  packets  successfully  transmitted  in  the  interval  ( 
T;,Ti+ 1  ]  ,  these  are  the  packets  that  arrived  during  the  interval  [  T,  -  l,T)+i  -  1  ).  The 
sequence  /?,  ;  i  >  1,  is  a  renewal  process  since  Qi ;  i 1  1,  is  a  sequence  of  nonnegative 
i.i.d  random  variables.  Furthermore,  the  delay  process  Dn  ;  n  >  1,  is  regenerative  with 

respect  to  the  renewal  process  Rj ;  i  >  l,  with  regeneration  cycle  Q  j .  From  the  regenei- 

Q\ 

tive  theorem  [20],  we  conclude  that  if  Q  =  E(Q\)  <  00  and  W  =  E  {  £  DJ  <  «,  then 

i  =  i 

there  exists  a  real  number  D  such  that, 

D  =  limn-1  =  lim  (££>;)  =  (HI.3) 

n~*°°  i=i  L? 

The  first  and  second  equalities  in  (3)  above,  are  w.p.  1  .  In  addition,  since 
P  (Qi  =  1)  >  0,  the  distribution  of  Q\  is  aperiodic  and  there  exists  a  random  variable  D„ 
such  that  the  sequence  Dn  ;  n  >  1  converges  in  distribution  to  D„.  £><«,  represents  the 
steady  state  delay  induced  by  the  algorithm  and  its  mean  satisfies  the  equality 

£(£>„)  =  -jj  (ffl.4) 

The  quantity  D  will  be  referred  to  as  the  mean  packet  delay.  From  (4)  we  observe  that  the 
mean  packet  delay  can  be  determined  by  computing  the  quantities  of  the  right  hand  side 
of  the  equality.  In  Appendix  B  we  develop  two  systems  of  linear  equations  whose  solu¬ 
tion  may  be  used  to  compute  the  mean  cycle  length  Q  and  the  mean  cumulative  delay  W. 

In  Table  III.l,  we  include  the  computed  upper  and  lower  bounds,  D“  and  D;  respec- 
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proposed  and  the  Capetanakis's  dynamic  algorithms. 


Proposed 

algorithm 

Capetanakis  dynamic 
algorithm 

X 

Dl 

Du 

Dl 

Du 

0.02 

1.562 

1.563 

1.563 

1.564 

0.06 

1.708 

1.716 

1.713 

1.719 

0.10 

1.888 

1.917 

1.903 

1.921 

0.16 

2.257 

2.363 

2.308 

2.362 

0.20 

2.607 

2.812 

2.712 

2.809 

0.24 

3.103 

3.467 

3.308 

3.476 

0.30 

4.412 

5.197 

4.976 

5.365 

0.32 

5.162 

6.170 

5.973 

6.501 

0.36 

7.941 

9.665 

9.798 

10.883 

0.38 

11.008 

13.398 

14.121 

15.855 

0.40 

18.262 

22.024 

24.427 

27.736 

0.42 

57.354 

67.665 

78.530 

90.212 

Table  III.l 


Upper  and  Lower  Bounds  on  Steady-State  Expected  Delays 


Regarding  the  throughput  X*  and  the  optimal  window  size  A*,  the  following  results  have 
been  found. 


Proposed  Algorithm: 

X  =  0.4295 

A'  =  2.33 

Capetanakis’  Dynamic  Algorithm: 

X *  =  0.4295 

A'  =  2.677 

Table  ni.2 


Throughputs  and  Optimal  Window  Sizes 


» 

From  Table  III. 2,  we  observe  that,  under  the  Poisson  user  model,  the  algorithm  in 
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algorithm,  but  uses  a  smaller  window  size.  From  Table  m.l,  we  observe  that  the  two 
algorithms  induce  practically  identical  delays  for  Poisson  rates  in  (0,  0.30),  while  for 
Poisson  rates  in  (0.30, 0.42],  the  proposed  algorithm  induces  lower  delays. 


Remarks  It  may  seem  surprising  that  the  algorithm  in  this  chapter  attains  the  same 
throughput  with  that  of  the  Capetanakis’s  dynamic  algorithm,  since  the  expected  lengths 
Lk  in  (2a)  are  bounded  by  a  quadratic  function  of  k,  while  the  corresponding  lengths  for 
the  Capetanakis’s  algorithm,  Lk  ,  are  bounded  by  linear  a  function  of  k.  An  intuitive 
explanation  is  as  follows:  The  proposed  algorithm  is  on  the  average  faster  than  the 
Capetanakis’s  algorithm,  when  the  collision  multiplicity  is  2,  since  Li  -  4.5  and  Li  =  5. 
For  collision  multiplicities  higher  than  2,  the  proposed  algorithm  is  on  the  average  slower 
than  Capetanakis’s  algorithm  since  Lk>  Lk  ,  k>  3.  Furthermore,  the  probability  of  a 
higher  than  two  collision  multiplicity  for  the  proposed  algorithm,  is  smaller  as  compare^ 
to  the  corresponding  probability  for  the  Capetanakis’s  algorithm  (this  is  because  the  pro¬ 
posed  algorithm  uses  a  smaller  optimal  window).  The  advantage  of  the  proposed  algo¬ 
rithm  for  collision  multiplicities  equal  to  2,  is  balanced  by  the  advantage  of  the 
Capetanakis’s  algorithm  for  collision  multiplicities  greater  than  2. 


HI.3  PERFORMANCE  UNDER 
FEEDBACK  ERRORS  AND  OPERATION  UNDER  LIMITED 
FEEDBACK  SENSING 

In  this  section,  we  study  two  important  characteristics  of  the  algorithm.  Its  perfor¬ 
mance  under  feedback  channel  errors,  and  its  operation  and  under  limited  feedback  sens¬ 
ing. 
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III.3.1  Performance  Under  Feedback  Errors 

Let  us  assume  that  due  to  noisy  conditions,  the  following  types  of  feedback  errors 
may  occur:  With  probability  £  an  empty  slot  may  be  perceived  by  the  users  as  a  collision 
slot.  Also,  with  probability  5  a  slot  occupied  by  a  single  transmission  may  be  perceived 
by  the  users  as  a  collision  slot.  We  assume  that  a  collision  slot  is  always  recognized 
correctly  by  the  users.  We  consider  the  case  where  the  probabilities  £  and  8,  are  known  a 
priori.  Then  given  E  and  8,  the  window  size  A  is  optimized  for  throughput  maximization. 
We  performed  throughput  analysis,  (the  details  are  included  in  Appendix  B),  for  both  the 
proposed  algorithm  and  the  Capetanakis’s  dynamic  algorithm,  [11].  We  exhibit  the 
results  in  Table  III.3.  From  Table  III.3,  we  conclude  that  the  proposed  algorithm  is  very 
insensitive  to  feedback  errors.  Even  for  the  practically  extreme  case  £  =  8  =  0.1,  the 
throughput  is  almost  90%  of  its  value  in  the  error  free  case.  We  also  conclude  that  the 
proposed  algorithm  allows  operation  (positive  throughput)  as  long  as  £  <  1  and  8  <  1, 
while  if  £>0.5  the  throughput  for  the  Capetanakis’s  algorithm  is  zero.  We  notice  for 
example  that  for  £  =  0.5  and  8  =  0,  the  proposed  algorithm  attains  throughput  equal  to 
A.*  =0.325. 


e 

5 

X  proposed  alg. 

r  Cap. 

0.00 

0.00 

t  0.4295 

0.4295 

0.00 

0.01 

0.4248 

0.4258 

0.00 

0.10 

0.3873 

0.3920 

0.00 

0.20 

0.3463 

0.3535 

0.00 

0.40 

0.2655 

0.2731 

0.00 

0.50 

0.2251 

0.2310 

0.01 

0.00 

0.4272 

0.4272 

0.10 

0.00 

0.4117 

0.4043 

0.20 

0.00 

0.3930 

0.3706 

0.40 

0.00 

0.3503 

0.2329 

0.45 

0.00 

0.3382 

0.1524 

0.50 

0.00 

0.3250 

0.0000 

0.60 

0.00 

0.3050 

0.0000 

0.70 

0.00 

0.2750 

0.0000 

0.80 

0.00 

0.2280 

0.0000 

0.90 

0.00 

0.1700 

0.0000 

0.10 

0.10 

0.3706 

0.3672 

0.20 

0.20 

0.3139 

0.2972 

0.30 

0.30 

0.2589 

0.2166 

0.40 

0.40 

0.2064 

0.1205 

0.30 

0.50 

0.1885 

0.1511 

0.30 

0.70 

0.1183 

0.t  "86 

0.90 

0.90 

0.0105 

0.0000 

Table  III .3 

Throughputs  as  a  function  of  e  and  5.  Window  sizes  optimized  for  every  pair  (e  ,  8) 

In  some  systems,  the  probabilities  e  and  5  may  not  be  known  a  priori.  In  this  case, 
we  designed  the  algorithm  subject  to  the  assumption  of  error  free  feedback.  The 
corresponding  results  are  shown  in  Table  III.4. 
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e 

5 

X  prop.  aig. 

X'  Cap. 

0.00 

0.00 

0.429 

0.00 

0.01 

0.424 

0.00 

0.10 

0.386 

0.00 

0.40 

0.255 

0.262 

0.01 

0.00 

0.427 

0.427 

0.10 

0.00 

0.408 

0.401 

0.20 

0.00 

0.377 

0.355 

0.3 

0.00 

0.320 

0.248 

0.3305 

0.00 

0.289 

0.000 

0.37 

0.00 

0.214 

0.000 

0.3787 

0.00 

0.000 

0.000 

0.10 

0.10 

0.363 

0.361 

0.10 

0.20 

0.318 

0.318 

0.10 

0.30 

0.273 

0.272 

0.20 

0.10 

0.328 

0.311 

0.20 

0.20 

0.279 

0.266 

0.20 

0.30 

0.232 

0.221 

0.30 

0.10 

0.261 

0.192 

0.30 

0.20 

0.205 

0.138 

0.30 

0.30 

0.153 

0.088 

Table  III.4 

Throughputs  as  a  function  of  £  and  8.  Window  sizes  unchanged  for  every  pair  (E  ,  5) 

From  Table  IEL4,  we  observe  that  the  proposed  algorithm  is  better  than  the  Capetanakis’s 
algorithm.  For  example  the  maximum  e  value  for  which  the  proposed  algorithm  is  stable 
is  0.378,  while  the  Capetanakis’s  algorithm  becomes  unstable  for  e  >  0.33. 

in.3.2  Operations  Under  Limited  Feedback  Sensing 

Under  limited  sensing,  it  is  required  that  each  user  monitors  the  channel  feedback 
from  the  time  he  generates  a  packet,  until  the  time  his  packet  is  successfully  transmitted. 
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our  objective  is  to  prevent  transmission  of  new  arrivals  until  the  current  collision  resolu 
tion  process  has  been  completed.  The  latter  can  be  accomplished,  provided  that  each 
user  with  a  new  arrival  can  decide,  in  a  finite  number  of  slots,  whether  a  collision  resolu¬ 
tion  process  is  in  progress.  A  user  with  a  new  arrival,  who  observes  a  C  slot,  decides  to 
wait  because  he  deduces  that  collision  resolution  is  in  progress.  Two  consecutive  NC 
slots  correspond  to  either  two  consecutive  CRIs  of  length  one,  or,  to  the  end  of  a  CR1 
which  started  with  a  collision  slot.  Therefore,  if  two  consecutive  NC  slots  are  observed, 
the  user  decides  that  the  system  is  empty.  Summarizing,  under  limited  feedback  sensing 
the  algorithm  can  be  modified  as  follows: 

The  window  size  remains  the  same  as  in  the  case  of  continuous  feedback  sensing.  The 
window  slides  through  the  unexamined  interval  from  the  current  time  to  the  past  (since  a 
user  with  a  new  arrival  has  no  information  regarding  the  events  that  happened  prior  to  the 
arrival  time).  Its  edge  is  maintained  one  slot  before  the  current  time,  (see  also  Figure  111.1 
).  Within  a  CRI,  the  algorithm  operates  exactly  as  under  continuous  feedback  sensing. 

For  very  light  input  rate  a  user  with  a  new  arrival  will  observe  an  empty  slot  (arrival 
slot)  with  high  probability.  He  then  waits  for  the  next  slot  (observation  slot)  which  will 
be  also  empty  with  high  probability.  The  observation  slot  being  empty,  the  user 
transmits  in  the  subsequent  slot  (transmission  slot).  Based  on  the  light  input  rate  condi¬ 
tion,  we  silently  assume  that  no  collision  occurs  over  the  transmission  slot.  Because 
Poisson  arrivals  over  a  time  interval  are  uniformly  distributed,  we  may  infer  that  the  new 
packet  will  show  up  in  the  middle  of  the  arrival  slot;  i.e.,  delay  over  the  arrival  slot 
equals  0.5.  Adding  up  the  delays  over  the  arrival,  observation,  and  transmission  slots,  we 


41 


As  the  input  rate  increases,  a  user  with  a  new  arrival  observes  a  collision  slot  with 
high  probability,  and  waits  for  its  resolution  (occurrence  of  two  consecutive  NC  feed¬ 
backs).  Therefore,  under  heavy  traffic  conditions,  a  new  user  can  determine  the  status  of 
the  channel  immediately  upon  arrival  (without  the  waste  of  any  observation  slots).  Con¬ 
sequently,  the  total  delay  equals  that  under  continuous  feedback  sensing.  The  throughput 
of  the  algorithm  remains  identical  to  that  under  continuous  feedback  sensing  (since  the 
window  size,  and  the  collision  resolution  algorithm  are  the  same  in  both  cases). 
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U:  unresolved  Interval  ct:  current  time 

R:  resolved  Interval 


Figure  III.  1 

Window  selection  in  the  limited  sensing  case. 
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m.4  THE  OUTPUT  TRAFFIC  INTERDEPARTURE  DISTRIBUTION 

In  this  section,  we  analytically  evaluate  the  steady-state  distribution  of  the  distance 
between  two  consecutive,  successful  transmissions,  when  the  proposed  algorithm  is 
employed.  The  algorithm  generates  an  output  traffic  process  with  memory.  Conse¬ 
quently,  our  computations  correspond  to  the  First  order  distribution  from  this  process. 
The  first  order  distribution  together  with  a  memoryless  assumption,  may  be  used  as  an 
approximation  of  the  actual  output  traffic  process,  when  one  studies  interconnected  sys¬ 
tems  which  employ  the  algorithm  for  their  internal  transmissions.  Typical  example  is  the 
multi-hop  problem,  where  part  of  the  output  (internally  successfully  transmitted  packets) 
of  a  random  access  system  has  to  be  transmitted  through  another  random  access  system. 

The  methodology  utilizes  the  regenerative  character  of  the  output  traffic  process. 
We  define  the  sequence  {PJi>\  as  follows:  Each  P,  is  a  collision  resolution  point, 
(CRP),  which  follows  a  slot  containing  a  successful  transmission  and  at  which  the  lag 
equals  one.  P \  is  the  first  after  zero  such  CRP,  and  for  every  <>1,  P,+]  is  the  first  after 
Pt  such  CRP.  Let  S,f  i>l  denote  the  number  of  i  (  '.ssful  transmissions  in  (0,  P,],  and 
let  dn  denote  the  distance  between  the  (n-l)-th  and  the  n-th  successful  transmission. 
Then,  5,,  i>\  is  a  renewal  process,  and  the  process  d ^  n>  1  is  regenerative  with  respect  to 
it.  Let  us  define,  C,  =  Si+1 -5,,  i>l.  Then  C,  denotes  the  number  of  successful 
transmissions  in  the  interval  (P,,  P,+IJ,  where  this  interval  will  be  called  the  i-th  cycle. 
Let  us  define, 

f  1  ,  if  dn  =  s 

,  otherwise  (III. 5) 

H  =  E{Pl+]-PJ  (III. 6) 

From  the  regenerative  theorem  (201,  we  conclude  that  if  C  =  E (C \}<  »,  then. 
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lim  N'1  £  {n(s)=  lim  N~l  E{%  In{s)}  =  C~x  E{£  In(s)}  (ffl.7) 

N^°°  n=l  N~*oe  n=l  n=l 

jwhere  the  first  equality  in  (7)  is  w.p.  1.  We  also  have  that, 

C  =  XH  (HI. 8) 

;  where  A.  denotes  the  intensity  of  the  Poisson  input  traffic.  In  addition,  since 
P(C  i  =  1)  >  0,  the  distribution  of  C  i  is  aperiodic  and  there  exists  a  random  variable  dx, 

such  that  the  sequence  dn,  «= 1,2,...  converges  in  distribution  to  d^.  Then,  dx 

represents  the  steady  state  interdeparture  distance  induced  by  the  algorithm,  and  its  distri¬ 
bution  satisfies  the  equality, 


,  C| 

P{d„=s)  =  C~x  In(s)}  (DI.9) 

n  =  1 

C, 

The  finiteness  and  the  computation  of  the  quantities  C  and  E{  £  in  (9)  are  related 

n  =  l 

to  the  existence  and  computation  of  appropriate  solutions  to  infinite-dimensionality 
linear  systems.  The  details  of  the  analysis  are  presented  in  Appendix  B.  In  Table  HI. 5,  we 
include  the  computed  upper  and  lower  bounds  on  the  probability  P(d„  =  s).  In  Figure 
III. 2,  we  plot  the  lower  bounds  against  s,  for  various  Poisson  input  rates  X. 
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X  =0.1 

X  =0.4 

s 

pl 

r  s 

Pus 

P‘s 

pu 

1  s 

1 

0.1420 

0.1427 

0.4702 

0.4728 

2 

0.0816 

0.0832 

0.2000 

0.2048 

3 

0.0704 

0.0739 

0.0998 

0.10806 

4 

0.0641 

0.0696 

0.0612 

0.0787 

5 

0.0537 

0.0603 

0.0401 

0.0603 

6 

0.0502 

0.0591 

0.0280 

0.0397 

7 

0.0420 

0.0503 

0.0196 

0.0264 

8 

0.0364 

0.0452 

0.0099 

0.0173 

9 

0.0332 

0.0431 

0.0057 

0.0094 

10 

0.0265 

0.0393 

0.0013 

0.0062 

Table  III.5 


Upper  and  lower  bounds  on  the  interdeparture  distribution  for  Poisson  rates 

X  =0.1  and  X  =0.4 


From  Table  113.5  and  Figure  III. 2,  we  conclude  the  following: 

(1)  For  small  X  values  (X  <  0.1),  the  interdeparture  distribution  is  close  to  the  Bernoulli 

distribution  with  parameter  p  =  Xe~ \  In  particular,  denoting  by  Ps^.P(dao=s),  we 

have  Ps=  p  (1  -p)s~l ,  for  s  >  2.  The  probability  P  i  however,  is  larger  than  the  Bernoulli 
parameter  p.  An  intuitive  explanation  of  the  latter  is  as  follows:  For  small  X  values,  sin¬ 
gle  arrivals  in  two  consecutive  slots  occur  with  probability  p2  =  (Xe~x)2  ~  X2.  TTie  pro¬ 
bability  of  a  collision  slot  is  approximately  equal  to  2~l  X2  e~x  ~  2~x  X2.  Therefore, 
under  light  traffic  single  arrivals  in  two  consecutive  slots  contribute  two  thirds  of  P  \ , 
^  while  the  remaining  one  third  is  due  to  consecutive  departures  at  the  end  of  a  collision 

resolution  interval. 

(2)  As  the  rate  of  the  Poisson  input  traffic  increases,  the  interdeparture  distribution 

I  induced  by  the  algorithm  deviates  further  from  the  Bernoulli  distribution.  In  that  case  the 

mass  of  the  distribution  accumulates  at  relatively  small  s  values  (see  Table  III. 5,  for 
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10 

X  =  0.4,  we  have  P  i  =  0.471  and  £  Ps  -  1. 

5  =  1 


Remark  Our  results  show  that  it  is  generally  wrong  to  conjecture  Bernoulli  interdeparture 
distribution. 
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Figure  III.2 
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Lower  bounds  on  the  interdeparture  distribution. 
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III.5  CONCLUSIONS 

We  presented  a  simple  window  random  access  algorithm  for  systems  with  binary, 
collision  versus  noncollision,  feedback.  We  analyzed  the  algorithm  assuming  the  Pois¬ 
son  user  model,  and  continuous  feedback  sensing.  In  addition  to  the  throughput  and 
delay  analyses,  we  studied  the  effect  of  feedback  errors  on  the  throughput  of  the  algo¬ 
rithm.  We  also  evaluated  the  output  traffic  interdeparture  distribution.  Both  the  pro¬ 
posed  and  the  Capetanakis’s  algorithm  attain  the  same  throughput.  In  terms  of  induced 
delays  and  feedback  error  insensitivity,  the  proposed  algorithm  is  better  than  the 
Capetanakis’s  dynamic  algorithm.  Furthermore,  the  proposed  algorithm  can  be  easily 
modified  to  operate  under  limited  feedback  sensing,  and  its  output  traffic  interdeparture 
distribution  can  be  analytically  evaluated. 
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CHAPTER  IV 


DELAY  DISTRIBUTION  ANALYSIS  OF  WINDOW  RANDOM-ACCESS  ALGORITHMS 


IV.  I  INTRODUCTION 

Window  Random-Access  algorithms  constitute  an  important  class  of  Multiple- 
Access  algorithms;  they  are  distributive  and  attain  high  throughputs  and  low  delays  by 
controlling  the  number  of  simultaneously  transmitting  users.  The  throughput  analysis  of 
algorithms  in  this  class,  is  a  relatively  easy  task.  The  delay  analysis,  however,  presents 
difficulties,  due  mainly  to  the  variable  window  sizes  and  the  complicated  state  space  that 
some  of  these  algorithms  create.  These  restrictions  prohibit  the  application  of  results 
from  standard  queueing  theory  in  the  delay  analysis. 

Many  attempts  for  the  delay  analysis  have  been  made.  In  [23],  the  class  of  algo¬ 
rithms  with  constant  window  size  was  considered,  and  upper  bounds  on  the  expected 
delays  were  developed.  In  [24],  [20],  methods  for  the  computation  of  bounds  on  the 
moments  of  the  delays  were  presented.  A  method  for  the  computation  of  delay  distribu¬ 
tion  for  constant  window  size  algorithms  appears  in  [25].  The  method  in  [25]  relies  on  a 
clever  decomposition  of  the  delay  process,  which  allows  the  application  of  results  from 
standard  queueing  theory.  The  computation  of  the  delay  distribution  for  variable  window 
size  algorithms,  however,  remains  an  open  problem.  One  possible  approach  is  to  com¬ 
pute  bounds  on  the  moments  of  the  delays  as  in  [24]  or  [20],  which  can  then  be  used  for 
an  approximate  evaluation  of  the  delay  distribution.  This  approach,  however,  is  not  com- 

. nro/'tio.al 
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In  this  chapter,  we  show  that  the  methodology  employed  in  [20],  can  be  extended  to 
provide  bounds  on  the  distribution  of  the  delays.  The  quantities  of  interest  are  related  to 
the  solution  of  a  denumerable  system  of  linear  equations.  Methods  for  the  computation 
of  the  constant  terms  and  the  coefficients  of  the  unknowns  of  the  system  are  developed. 
The  methodology  is  applied  to  the  delay  distribution  analysis  of  both  the  Capetanakis 
Window-Access  algorithm  with  binary  feedback  and  the  Part-and-try  algorithm  with 
binary  feedback.  It  can  also  be  applied  directly  to  other  Window  Random-Access  algo¬ 
rithms  with  different  feedback.  An  interesting  result  of  the  analysis  is  that  as  the  arrival 
rate  increases,  the  tails  of  the  distribution  become  longer,  but  the  median  grows  much 
slower  than  the  expected  delay. 

IV.2  MODEL  SPECIFICATION 

We  consider  a  single  slotted  channel  that  is  being  accessed  by  a  number  of  indepen¬ 
dent  packet  transmitting  users.  The  length  of  a  packet  is  equal  to  the  length  of  a  slot,  and 
packet  transmission  may  start  only  at  the  beginning  of  a  slot.  Simultaneous  transmission 
of  more  than  one  packets  in  the  same  slot,  results  in  complete  loss  of  the  information 
included  in  the  involved  packets.  The  latter  event  is  referred  to  as  a  "collision”  event.  At 
the  end  of  each  slot,  all  users  receive  a  feedback  that  provides  some  information  about 
the  channel  activity  in  that  slot.  Common  types  of  feedback  are  the  binary  C-NC  (colli¬ 
sion  versus  noncollision)  feedback,  and  the  ternary  0-1 -C  (empty  versus  success  versus 
collision)  feedback.  To  resolve  the  collision,  the  users  follow  the  rules  of  a  Random- 
Access  algorithm.  The  algorithm  is  implemented  by  each  user  in  a  distributed  fashion, 
using  only  the  available  feedback.  The  cumulative  packet  generating  process  is  assumed 
to  be  Poisson  with  rate  X  packets  per  slot. 
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We  assume  that  a  Window  Random-Access  algorithm  is  employed,  whose  basic 
operating  characteristics  are  the  following  (see  Figure  IV.  1):  Suppose  that  at  the  begin¬ 
ning  of  slot  v  all  packet  that  arrived  before  time  rv<v  have  been  successfully  transmitted, 
and  there  is  no  information  concerning  the  packets  that  may  have  arrived  in  the  interval 
[fv,v),  (i.e.,  the  distribution  of  the  interarrival  times  of  the  packets  in  [rv,  v)  is  the  same  as 
the  one  assumed  originally).  The  beginning  of  such  a  slot  v  is  called  a  Conflict  Resolu¬ 
tion  Point  (CRP).  The  time  difference  dv  =v-tv  is  referred  to  as  the  lag  at  v.  In  slot  v, 
the  users  that  generated  packets  in  the  interval  [rv,  rv+xv),  where  Xv  =  min  (dv,A),  are 
allowed  to  transmit;  A  is  a  parameter  to  be  properly  chosen  for  throughput  maximization. 
After  a  random  number  of  slots  and  following  the  rules  of  the  algorithm,  another  CRP, 
v  ,  is  reached,  with  a  corresponding  tv  >  tv.  All  the  packets  that  have  been  generated  in 
the  interval  [rv  rv  ),  have  been  successfully  transmitted  in  the  interval  [v,v  ).  The  inter¬ 
vals  [v,v  ),  [rv,  tv+zv),  [rv,  fv  )  are  called  conflict  resolution  interval,  transmitted  inter¬ 
val,  and  resolved  interval,  respectively.  The  length  of  xv,  is  called  the  window  size  at  time 
v  .  Clearly,  the  window  size  varies  with  time,  and  its  maximum  size  is  A.  Note  also,  that 
the  length  of  the  conflict  resolution  interval  is  one,  if  and  only  if  there  are  at  most  one 
packets  in  the  transmitted  interval. 

Algorithms  that  operate  as  described  above,  are  the  Capetanakis  Window  Random- 
Access  algorithm  [11]  and  the  Part-and-Try  algorithm,  under  either  binary  C-NC  feed¬ 
back  [29],  or,  ternary  0-1 -C  feedback  [11], 

IV.3  STEADY  STATE  DELAY  DISTRIBUTION  ANALYSIS 

Let  packets  be  labeled  1,2,3,...  according  to  the  order  of  their  arrival  instants.  The 
delay  Dn  experienced  by  the  n-th  packet  is  defined  as  the  time  difference  between  its 


t 


:  generation  instant  of  the  j-th  packet 


\ 


:  successful  transmission  of  the  j-th  packet 


Figure  IV.  1 

Illustration  of  relationships  among  certain  random 
variables  related  to  the  operation  of  the  algorithm. 
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steady  state  distribution  of  Dn,  when  it  exists. 

Let  v,  ;  i>l  be  the  sequence  of  successive  CRPs  and  let  d,-  be  the  lag  at  v,.  The 
sequence  d,  ;  i2l  is  a  Markov  chain  with  state  space  F.  For  most  of  the  existing  Window 
Random-Access  algorithms,  F  is  a  denumerable  subset  of  the  interval  [1,®°).  Let  T j=l, 
dj  =  l,  and  define  7,+1,  as  the  first  CRP  after  7,,  at  which  the  length  of  the  lag  equals  one. 
From  the  description  of  the  algorithm  it  can  be  seen  that  the  induced  delay  process  proba¬ 
bilistically  restarts  itself  at  the  beginning  of  each  slot  7/,  /= 1,2,...  The  interval  [7,-,7,+1) 
will  be  referred  to  as  the  i-th  session.  Note  that  the  sessions  have  lengths  that  are  i.i.d. 
random  variables. 

Let  /?,,  t =1,2,...,  denote  the  number  of  packets  successfully  transmitted  in  the  inter¬ 
val  (0,7J;  (note  that  R,  also  represents  the  number  of  packets  arrived  during  the  interval 
[0,  7,-1),  since  7,  is  a  CRP  at  which  the  lag  is  one).  Then,  C,=Ri+i~R />1,  is  the 
number  of  packets  successfully  transmitted  in  the  interval  (7,7,+i  ]  -  these  are  the  pack¬ 
ets  that  arrived  during  the  interval  [7,-1,  T;+i-l).  The  sequence  j>1,  is  a  renewal 
process,  since  C,,  t>l,  is  a  sequence  of  nonnegative  i.i.d.  random  variables.  Further¬ 
more,  the  delay  process  D„,  n>l,  is  regenerative  with  respect  to  the  renewal  process 
/?,,  />1,  with  regeneration  cycle  Cj. 

Let 


,  ,J  1  if  Dn<s 

n  5  |  0  otherwise 


From  the  regenerative  theorem  [20],  we  conclude  that  if  C  =E(C\ )  <  then 

c, 

E(ZIn(s)) 

lim  -Jr  Z  In(s)  =  Hm  77F ( I  In(s))  =  — — - 


(IV.  1) 


N->~  /V 


n=\ 


n  =  l 


distribution  of  C\  is  aperiodic  and  there  exists  a  proper  random  variable  D^,  such  that  the 
sequence  Dn\  n=l,2,...  converges  in  distribution  to  D„.  D„  represents  the  steady  state 
delay  induced  by  the  algorithm  and  its  distribution  satisfies  the  equality 


c, 

£(£/„(*)) 

P(D„<s)  =  — ^ - 


(IV.2) 


From  (2)  we  observe  that  the  steady  state  distribution  of  the  delays  can  be  determined  by 

computing  the  quantities  of  the  right  hand  side  of  the  equality.  In  [20]  it  was  shown  that 

the  finiteness  and  the  computation  of  C  is  related  to  the  existence  and  the  computation  of 

an  appropriate  solution  to  an  infinite  system  of  linear  equations.  In  this  section  we  will 

C, 

show  that  the  same  is  true  for  the  quantity  E(  £/„(s) ). 

n=l 

The  following  definitions  will  be  used  in  the  sequel. 

/  :  Length  of  a  conflict  resolution  interval 

5  :  Length  of  a  resolved  interval 


x  :  Window  size. 


E  ( X  lx)  :  Expected  value  of  the  random  variable  X,  given  that  the  window  size  is 

x. 


p(x,r/x)  :  The  probability  that  the  conflict  resolution  interval  has  length  x  and  the 
resolved  interval  has  length  r,  given  that  the  window  size  is  x. 

p(x/x)  :  The  probability  that  the  conflict  resolution  interval  has  length  x,  given 

that  the  window  length  is  x. 

hj  :  Number  of  slots  needed  to  reach  a  CRP  with  lag  1  given  that  we  start 
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kd(s)  :  Number  of  successfully  transmitted  packets  with  delay  less  than  s,  in  the 

interval  hd. 

mXid(s)  :  Number  of  successfully  transmitted  packets  with  delay  less  than  s  during 
a  conflict  resolution  interval,  given  that  the  window  size  is  t  and  the  lag  is 
d. 

nx(s)  :  Number  of  successfully  transmitted  packets  with  delay  less  than  s,  dur¬ 

ing  a  conflict  resolution  interval,  given  that  the  window  size  is  t  and  the 
lag  is  x.  That  is,  nx(s)=mxx(s). 

Let  us  also  define. 


Kd(s)  =  E(kd  (s)) 

M  x,d(s)  =  E  (mX'd(s)) 
Nx(s)  =  E(nx(s)) 
Hd  =  E{hd ) 

Note  that  by  definition, 

c, 

n  =  l 

Also, 


C  =kH  j 

Therefore,  the  determination  of  ^i(j)  and  H\,  will  permit  the  computation  of  the 
steady-state  distribution  of  the  delays. 

Consider  the  arrangement  of  Fig.  IV.  1.  The  delay  D  of  the  successfully  transmitted 
packet  1,  can  be  decomposed  as  follows: 

D  —  0 1  — X+02 

Therefore, 
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D  £  s  iff  0 1  +02  —  J* — ' d  +z 


But  0i  4-02  is  statistically  identical  to  the  delay  that  the  successfully  transmitted  packet 
experiences  if  the  transmitted  interval  is  T  and  the  lag  is  t.  The  last  observation  shows 
that  mxd(s)  is  identically  distributed  with  nx(s-d+x).  Observe  now  that 


kd(s)  = 


mx yd(s)  if  dv'  =  l 

mx  d(s)+kj  (s)  if  dy,=d  * 1 


and  that 


dv‘  —  dv — &+/,  x  — 


dy,  if  dv< A 
A  ifdu>A 


From  (4)  and  (5)  we  conclude  that 


(IV.4) 


(IV.5) 


Kd(s)  =  Md d(s)  +  '£Kd-r+x(s)P(x,r/d)  if  l^f<A,  deF 


(IV. 6a) 


Kd(s)  =  MA.d(s)  +  Yfid-r+x(s)p(x,r/A)  if  d>A,  deF 


(IV. 6b) 


Since  mx  d(s)  is  identically  distributed  with  nx(s-d+X),  equations  (6a)  and  (6b)  become, 


Kd(s)  =  Nd(s)  +  '£Kd-r+x(s)p(x,r/d)  if  l<d<A,  deF 


(IV. 7a) 


Kd(s)  =  N zis-d+A)  +  ^Kd-r+x^)p(x,rlA)  if  d> A,  deF  (iv.7b) 

ryx 

Equations  (7a),  (7b)  comprise  a  denumerable  system  of  linear  equations.  Of  interest  to  us 
is  the  element  K\(s)  of  a  particular  solution  of  this  system.  The  methodology  developed 
in  [20]  can  be  used  for  the  study  of  system  (7).  Note  that  the  coefficients  of  the 
unknowns  are  independent  of  s.  This  observation  represents  a  computational  advantage 
when  the  solution  to  (7)  is  approximated  by  the  solution  of  appropriate  finite  linear  sys¬ 
tem  of  equations,  [20],  In  this  case,  the  approximate  solution  can  be  represented  in  the 
form. 


55 


be  computed  once,  and  then  used  for  the  computation  of  the  approximate  solution  for 
various  values  of  s. 


We  now  proceed  in  the  development  of  an  initial  upper  bound  on  the  solution  of 
system  (7).  Following  the  methodology  in  [20],  such  a  bound  will  be  the  sequence 
Kd(s)  =  if  yu(s)>  CuCO,  can  be  determined  so  that  the  following  inequali¬ 

ties  are  satisfied 


Kd(s)^  Nd(s)  +  '£Kd-r+x(s)p(x,r/d)  =  Kd(s)  if  l<d<A,  deF 

r.x  (IV. 8a) 

x»l 

K°d(s)>NA(s-d+&)  +  ZK°d_r+x(s)p(x,r/A)  =  Kld(s)  if  d>A,  deF  (IV  gb) 

Substituting  Kd{s)  in  the  right  hand  side  of  inequalities  (8),  it  can  be  seen  that  if  deF, 

=  K°d(s)  +  Nd(s)  +  yu(s)(E(l/d)  -  E (8/d)  -  (1  +Xd)e~™)  -  Us)d+U)e-u 

if  1  <d  <A  (IV. 9a) 

Kd{s)  =  K%)  +  N  A{s-d+b)-yu{s){E(§i&)-E{l/&))  ifd>A  (IV.9b) 

Observe  now  that  N A(s)  is  an  increasing  function  of  s.  Therefore,  from  (9b)  we  conclude 
that, 


Kld{s)<K°d{s)  +  NA{s)-yu(s)(E(b/&)-E(l/A))  ifd>A  (IV.  10) 

From  (10)  we  conclude  that  if  E{V A)  <  £(8/A),  the  condition  for  stability  of  the  system, 
inequalities  (9b)  are  satisfied  if 


y u(s)= 


Na(s) 

E  (5/A>-£  (//A) 


(1V.11) 


With  this  value  of  y u(s),  it  can  be  seen  that  inequalities  (9a)  are  satisfied  if 
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V(d)  = 


Nd(s )  +  yu(s)[E(Ud)  -  £(5/d)  -(l+Xd)e~u] 
(1  +kd)e~u 


(IV.  13) 


From  the  above  discussion  we  conclude  that  the  solution  to  system  (7)  satisfies  the  ine¬ 
qualities 


Kd(s)<yu(s)d  +  t,u(s),  deF 


(IV.  14a) 


where  yu(s),  £ u(s ),  are  given  by  equations  (1 1),  (12)  respectively.  The  uniqueness  of  the 
solution  is  guaranteed  by  the  same  techniques  as  in  [20].  If  we  use  a  similar  method  for 
the  development  of  a  lower  bound,  we  find  that 


Y i(s)d  +  ^,(s)  =  Kd(s),dsF  (TV.  14b) 

where 


Yi(s)=0  and  C/(J)=  inf  {(Nd(s)/((l+kd)e-u)} 

l<d<A 


Bounds  on  H  x  are  given  in  [20],  formulas  (18),  (19). 

As  is  explained  in  [20],  the  bounds  (14)  can  be  used  to  further  improve  the  bounds 
on  ATi(s).  To  proceed  further,  however,  the  computation  of  the  quantities 
E(lld),  E(5/d),  p{x,r/d),  and Nd(s)  is  necessary.  In  section  4,  we  present  a  method  for 
the  computation  of  Nd(s)  for  the  Capetanakis  dynamic  algorithm,  and  then  proceed  in  the 
computation  of  tight  upper  and  lower  bounds  on  the  quantities  of  interest.  In  section  5, 
we  present  a  method  for  the  computation  of  the  quantities  Nd(s)  andp (x,rld)  for  the 
Part-and-Try  algorithm,  and  then  we  develop  bounds  on  the  distribution  of  the  delays. 
Due  to  the  complicated  state  space  of  the  latter  algorithm,  the  development  of  tight 
bounds  for  high  input  rates  becomes  computationally  cumbersome. 
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IV.4  THE  CAPETANAKIS  ALGORITHM  WITH  C-NC  FEEDBACK 

A  complete  description  of  the  algorithmic  rules  can  be  found  in  [1 1],  In  this  algo¬ 
rithm,  the  resolved  interval  is  always  equal  to  the  window  size  i.e.,  x  =  5.  This  results  in 
a  significant  simplification  of  the  state  space  F,  when  the  maximum  window  size  A  is  a 
rational  number.  The  restriction  of  A  to  the  rational  numbers  facilitates  the  development 
of  tight  bounds  for  the  delay  distribution  and  does  not  represent  any  disadvantage  in 
practice.  Moreover,  if  m  packets  are  involved  in  a  conflict,  the  conflict  resolution  process 
depends  only  on  m  and  not  on  the  window  length  or  the  generation  time  of  each  of  the 
packets.  The  last  property  facilitates  the  development  of  efficient  methods  for  the  com¬ 
putation  of  the  quantities  of  interest. 

Since  x  =  5,  we  have  that,  £(8/d)  =  d  and 

0;l^se 

Formulas  for  the  computation  of  E  (lid),  p(x/d),  can  be  developed  by  following  the  rea¬ 
soning  in  [11],  In  the  next  section  we  develop  a  method  by  which  the  quantities  Nd(s) 
can  be  computed. 

IV.4. 1  The  Computation  of  N^(s). 

Let  M  be  the  number  of  packets  in  the  window  d  (see  Fig.  IV. 2).  Then, 

Nd(s)  =  E(nd(s))  =  E(E(nd(s)l  M))=  £  E (nd(s)  /  M  =  m)e~u (iv  15) 

m  =  l 

Let 


/  _  J  1  the  delay  of  the  i-th  packet  in  d  is  at  most  s 

i.mt  •'~|0  otherwise 


(IV. 16) 


r 

i 


58 


Figure  IV. 2 

Illustration  of  the  delay’s  decomposition  for 
the  Capetanakis’s  algorithm. 
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The  "i-th  packet"  in  (16),  is  the  i-th  packet  in  a  random  enumeration  of  the  m  packets 
(randomly  chosen  packet).  Then 

m 

nd(s)=ZAm(s)  (IV.  17) 

i=l 

Given  M,  the  generation  time  of  each  of  the  M  packets  is  uniformly  distributed  in  the 
window  d  and  independent  of  the  generation  time  of  the  rest  of  the  packets.  From  this 
observation  and  the  definition  in  (16)  we  conclude  that  the  random  variables 
J i,m(f  )>  1  <i£m,  are  identically  distributed  (although  not  independent).  Therefore,  from 
(17)  we  conclude  that 

E  ( nd(s )  /  M  =  m)  =  mE(JX  m(s)  /  M  =  m)  (IV. 18) 

Let  0i ,  (82),  be  the  delay  of  the  First  randomly  chosen  packet  before  (after)  the  initializa¬ 
tion  of  the  collision  resolution  process  i.e.  until  (after)  time  v  (Fig.  IV.2).  Since  9i  is  uni¬ 
formly  distributed  in  [0,  d),  we  have  that 

E(JLm(s)lM=m)=-jE(Jx.m(s)IM  =  m,Q:  =Q)dQ  (IV.19) 

do 

Since  the  conflict  resolution  process  is  independent  of  the  packet  generation  time  in  a 
window,  we  conclude  that  given  m,  the  random  variables  0lt  02,  are  independent.  There¬ 
fore. 

E(J i(s)  /  M  =  m,  0]  =  0)  =p  (01+0;  <  s  /  M  =m,Q\  =  9)  =  p(02  S  5-0  /  M  =  m)  (IV. 20) 
and 

,  * 

E(Jlm(s)l  M  =  m)  =  —  f  p  (02  <s-Q/  M  =  m)dQ  (IV.21) 

d  o 

Observe  now,  that  02  takes  only  positive  integer  values.  Let 

Pf)=p(Q1=q  !M  =  m)\  q= 1, 2,3,... 

then. 
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t*-®!  ,  , 

<s-Q/  M  =  m)  =  £  PT] 

q=\ 


(IV. 22) 


¥-» 

where,  [a]  denotes  the  integer  par*  of  a.  Substituting  (22)  in  (21),  we  conclude  that 


1  , 

E(Jhm(s)/M=m)=±j  £  P(qm)  dQ 

&  A  „  _  1 


(IV.23) 


u  <7  =  i 

Let  e=d-s+[s].  Then,  since  is  independent  of  0,  (23)  can  be  written  as  follows: 

,  Ii]  ,  ,  [J]-1  .  .  I*H<I  ,  ,  (■*]-[« l'l  ,  , 

£(/,,m(s)/A/=^)  =  ^((j-[5])S^m>+  I  P(”)+...+  £  P^m)+(e-[e])  £  ^m)) 

“  ij=l  </  =1  <7=1  <7=1 

From  (15),  (18)  and  (24),  we  observe  that  for  the  computation  of  Nd(s),  only  the  quanti¬ 
ties  Pq n)  are  needed.  In  Appendix  C,  we  provide  recursive  formulas  for  the  computation 
of  the  quantities  It  is  worth  noticing  that  formulas  (15),  (18)  and  (24),  arc  valid  for 
any  window  Random-Access  algorithm  that  has  the  properties  described  in  the  first  para¬ 
graph  of  section  3. 


IV.4.2  Development  of  bounds  on  K  i  (s),  H  i 


For  the  development  of  bounds  on  the  delays  we  chose  A=2.5  As  a  result,  the  state 
space  F  becomes, 

F  =  {  1,  1.5,  2.,  2.5 . } 

For  the  algorithm  considered  in  this  section,  the  maximum  throughput  is  achieved  for 
A=2.67  [11].  The  reduction  in  throughput  due  to  the  choice  A=2.5  is  not  significant  (less 
than  .1%). 

Since  5st,  equations  (7a),  (7b),  become: 


Kd(s)  =  Nd(s)  +  '£KAs)p(x/d)  ifd=l,  1.5,  2,  2.5 


(IV. 25a) 


Kd(s)  =  NA(s-d+A)  +  '£Kd^+x(s)p(x/A)  if  d=3,  3.5,  4,...  (iv.25b) 

X 

For  the  development  of  bounds  on  ^i(s)  we  followed  the  method  of  truncation  of  the 
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infinite  system  (25),  see  also  [20].  Specifically,  in  system  (25),  we  replaced  the 
unknowns  Kd(s),  ford>40  with  the  upper  bounds  in  (14a).  The  substitution  results  in  a 
finite  system  of  equations  whose  solution  is  an  upper  bound  to  the  solution  of  (25)  for 
d= 1,  1.5,  2,. ..,40.  For  the  development  of  lower  bounds,  we  replaced  the  unknowns 
Kd(s),  d> 40  with  the  lower  bound  in  (14b).  The  same  methodology  is  employed  for  the 
development  of  bounds  on  H\.  The  resulting  upper  and  lower  bounds  on  the  distribution 
of  the  delays  differ  by  .01  in  the  worst  case.  In  Figure  FV. 3,  we  provide  the  distribution 
of  the  delays  for  various  values  of  the  arrival  rate.  An  important  observation  is  that  as  the 
arrival  rate  increases,  the  tails  of  the  distribution  become  longer,  but  the  median  grows 
much  slower  than  the  expected  delays. 

IV.5  THE  PART-AND-TRY  ALGORITHM  WITH  C-NC  FEEDBACK 

A  detailed  description  of  the  algorithm  can  be  found  in  [29].  The  algorithm  has 
throughput  .45,  and  its  basic  difference  from  the  Part-and-Try  algorithm  with  ternary 
feedback  in  [11]  is  that  if  a  collision  is  followed  by  an  empty  slot,  the  packets  that  are 
involved  in  the  collision  retransmit  in  the  next  slot  (no  splitting  takes  place).  The  tech¬ 
niques  used  in  this  section,  can  be  easily  applied  to  the  analysis  of  the  algorithm  under 
ternary  feedback. 

IV.5.1  The  computation  of  NAs). 

Let  M  be  the  number  of  packets  in  the  window  d  (see  Fig  IV.4).  As  in  section  4.1, 
formula  (15)  holds.  Since  not  all  the  packets  in  a  window  are  successfully  transmitted, 
however,  we  need  to  modify  the  definition  o(Ji  m(s)  in  (16).  Let 
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(1  if  the  i  -th  randomly  chosen  packet  in  d  is  sucessfully 

transmitted  and  its  delay  is  at  most  s.  (IV.26) 

0  otherwise 

then,  as  in  (18), 

E  (n^is)  /  M  =m)  =  mE  {J\<m{s)  /  M  =m)  (IV. 27) 

and 

E(JXm(s)/M=m)=-!-  \  E(Jlm(s)l  M=m,  01=0V0  (IV.28) 

d  Jo 

Again,  02  takes  integer  values.  However,  and  02  are  not  independent.  To  proceed 
further,  we  define  the  events 

Ai  m(s)={  the  i-th  randomly  chosen  packet  is  successfully  transmitted  and  02<5.j 
aim(q)  =  { the  i-th  randomly  chosen  packet  is  successfully  transmitted  and  02=<7,  q  = 1,2, ...J 

Since  0  takes  only  integer  values,  we  have  that 

U) 

Ai,m(s)=^jaim(q)  (IV. 29) 

<?=i 

From  (28)  and  (29)  we  conclude  that 

E(JUm(s)IM=rn)  =  -\p(AUm(s-Q)l  M=m,  Qx=Q)dQ  (IV.30) 

do 

1  ^  [*-oi 

=  “jJ  L  P(ai.m(q)/ M=m,  0i=0)d0 

“o  <7=1 

Observe  now  that  p(ax ,m(q)  /  M=m,  0j  =0)  depends  on  0  through  the  ratio  0=0/d.  Let  us 
define 

f:(/>  m,ty)  =  p(aim(q)l  M=m,  Bx=$d)\  4>=-^-e[0,l)  (IV.31) 

d 

then,  the  following  equations  hold: 


t 


/t(0,m,<5))  =  0,  for  m  =  1,2 . and  4>e  [0, 1) 


(IV. 32a) 


/t(l,m,<j>)  =  0,  form=2,3,..,  and<>e[0, 1) 


/l  (?.*■,$)*•{ 

m£2 


j=° 


m- 1 

J  J 


2«m-l)) 


.5<<(><1 
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avj2b> 

(IV. 32c) 

(IV.  3  2d) 


Formulas  (32a),  (32b),  (32c),  are  clear.  Let  us  explain  formula  (32d). 

Assume  that  .5  <<{><  1.  Then,  the  packet  under  consideration  (packet  1  in  Figure  IV. 4), 
lies  in  the  left  hand  half  (l.h.h.)  of  the  window.  Let  02  be  the  delay  of  packet  1,  after  the 
first  collision.  Then, 


02  =  02+1 


The  location  of  packet  1  in  the  l.h.h.  of  the  window,  is. 


0i  =0!-d/2 


Therefore,  the  new  ratio  becomes 


,  0j 

0  =  TT  =  20-1 


(IV. 33) 


(IV.  34) 


d/2 


(IV.35) 


Let  Gj .m_i  be  the  set 


Gj  m ~\={j  of  the  m- 1  packets  (other  than  packet  1),  are  located  in  the  l.h.h.  of  the  window.} 


Note  that p(Gj%m.\)  = 


m- f 
L  J  J 


2  (m  1}.  From  (33),  (34)  and  (35),  we  derive  formula  (32d) 

for  .5  <<|>5  1  by  conditioning  on  G;m_ If  0<4><.5,  packet  1  is  located  in  the  right  hand 
half  (r.h.h.)  of  the  window.  In  this  case,  '*'e  note  that  if  in  the  l.h.h.  of  the  window  there 
are  more  than  one  packets,  then  packet  1  is  not  transmitted  during  the  collision  resolution 
process.  By  conditioning  again  on  the  G(m_i,  and  taking  into  account  the  last  observa- 
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Figure  IV.4 

Illustration  of  delay’s  decomposition  for  the  Part-and-Try  algorithm. 
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From  equations  (32)  we  conclude  by  induction  the  following  property  for  the  func¬ 
tion/^,™,  <j>): 

Property  1.  For  fixed  q  and  m,  the  function  f\(q,m,§)  is  simple  (i.e.,  it  takes  a  finite 
number  of  values),  and  left  continuous.  The  jumps  of  the  function  occur  at  the  points 

0X(q~X),  ....  ....  l-2“<«-l>. 

Taking  into  account  Property  1,  we  can  compute  E (J\%m(s)),  using  formula  (30). 

IV.5.2  The  Computation  of  p  ( x,r/d ). 

The  method  of  conditioning  on  the  number  of  packets  in  a  window  (applied  in  sec¬ 
tions  4.1  and  5.1),  does  not  seem  to  lead  to  easily  computable  recursive  formulas.  In  this 
section,  we  present  an  alternative  methodology  that  results  in  simple  recursive  formulas 
for  the  computation  of  p(x,r/d),  and  provides  insight  into  the  structure  of  these  probabili¬ 
ties. 

Observe  first  that  p(x,r/d)  depends  on  r  through  the  ratio  s=r/d.  Let 
f2(x,s,d)  =  p  (/= x,  5=sd,  M>2  /  d) 

Conditioning  on  the  events  (M=0J,  { M=l }  and  {M> 2),  and  observing  that  if  M=0,1  then 
/=1  and  S=d,  we  conclude  that 

p(x,r  /d)  =  5K(x-\)SK«s  /d)-\)e-{U)(l+\d)  +  f2(x,s,d)  (IV.36) 

where 

5*'(y )  =  |" 1  if>,={) 

K{y>  1  0  otherwise 

The  function  f2(x,s,d)  satisfies  the  following  recursive  formulas 
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I 


> 


fli.X-.Syd) 


f2(x-l,2s,d/2)  0<s<5 

'  f2(x-2,2s-\,dl2)e~<Ml/l)(\+'>ull2 )  -5<s<l 
+  8k(x -3)8*(s -l)e“(W)(Xd/2)2 


(IV. 37b) 


We  derive  equation  (37a)  for  0<Ly<  5  The  case  .5<s<l  can  be  derived  by  a  similar  reason¬ 
ing.  Let  M  be  the  number  of  packets  in  the  l.h.h  of  the  window.  Let  l  be  the  number  of 
slots  during  a  collision  resolution  process  after  the  first  collision.  Then, 

p  (l=x,  5=sd,  M22  /  d)=p  (l  +  l=x,  6=sd,  M>2 ,  M  >2/  d)  +  p  (l—x,  &=sd,  M>2,  M  <1  /  d) 

But  if  M  <1,  then  s>.5.  Therefore  p(l=x,  &=sd,  M>2,  M  <1  / d)  =  0  if  0.£sr<5.  It 
remains  to  observe  that 


p(l+  l=x,  fcsd,  M>2,  M'>2  I  d)=p  (/  =x-l,  8=2(sd/2),  M  >2  /  d) 
and  that  if  M  >2,  the  collision  resolution  process  after  the  first  collision  is  statistically 

identical  to  a  collision  resolution  process  that  started  with  a  window  of  size  d/2. 

From  (37a),  (37b)  we  conclude  by  induction  that  the  function  f2(XyS,d)  has  the  fol¬ 
lowing  property: 


Property  2  For  fixed  x>3,  and  for  any  d,  s  takes  a  finite  number  rx  of  values.  The 
sequence  rx  satisfies  the  following  recursions:  r3=r4=l,  rx=rx_i+rx_2\  x>5.  Let  be 
the  set  of  values  of  s  for  given  x.  Then,1 

Al={V,  A4={\/2},  Ax  =  (.5 A,_,  )U(.5^_2+.5) ;  t>5 

The  values  of  f2(x,s,d)  for  various  x  and  s,  can  be  computed  from  (37b).  The  probabili¬ 
ties  p(x,r/d),  are  easily  computed  from  (36). 
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IV.5.3.  The  Development  of  bounds  on  K  i 


The  state  space  F  for  the  Part-and-Try  algorithm  is  a  dense  subset  of  [1,  °°).  This 
property  complicates  the  development  of  finite  systems  of  linear  equations  whose  solu¬ 
tion  will  provide  upper  or  lower  bounds  to  the  quantities  K  \  (5),  H ! .  Taking  advantage  of 
the  structure  of  the  probabilities  p(x,r/d)  (see  section  5.2),  however,  we  can  proceed  as 
follows:  Let  B  be  a  finite  subset  of  F,  that  includes  the  state  d=l.  We  develop  the  follow¬ 
ing  finite  system  of  linear  equations: 


Yd  =  Nd(s)+  X  (Y uU)(d-r+x)  +  Zu(s))p(x,rld)+  £  Yd.r+Xp(x,r  Id) 

( d-r+x)€B ‘  ( d-r+x)eB 


if  1  <d  <  A,  deB 


(IV. 38a) 


Yd  =  N&(s-d+A)  +  £  (yu(s)(d-r+x)  +  1^(s))p(x,r/d)  +  £  Yd_r^xp(x,rl A) 

( d-r+x)€B*  (d—r+x)€B 


if  d  >  A,  deB 


(IV. 38b) 


yu(s)  and  ^u(s)  are  determined  from  (11),  (12).  Due  to  property  2  of  section  4.2,  it  is  sim¬ 
ple  to  find  for  a  given  deB,  the  values  of  x  and  r  such  that  ( d-r+x )  e  B.  The  summation 
over  the  infinite  set  Bc ,  can  be  computed  in  terms  of  E  {l/d),  E  (bid),  y u(s),  £„( s ),  and  the 
probabilities  p(x,rld)\  (d-r+x)  e  B.  The  solution  Yd;  deB  of  system  (38),  is  an  upper 
bound  to  the  solution  Kd(s)\  deB  of  system  (7),  see  also  [20],  Since  leB,  we  can  deter¬ 
mine  a  bound  on  ^i(r).  Similarly,  lower  bounds  on  Kj(s),  and  upper  and  lower  bounds 
on  H\  can  be  developed.  For  the  computations  we  used  the  set 
B  =  {  1,  1.125,...,  l+(/!:/8),...,  9 }  .  The  resulting  bounds  on  the  delay  distribution  are 
presented  in  Figure  IV. 5  for  \  =  .1,  .2,  .3.  For  higher  arrival  rates,  the  bounds  are  not 
tight  and  although  they  can  be  improved  by  enlarging  the  set  B,  the  computations 


become  cumbersome. 


Delay  distribution  for  the  Part-and -Try  algorithm 
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IV.6  CONCLUSIONS 

We  developed  a  method  for  the  computation  of  bounds  on  the  delay  distribution  of 
Window  Random-Access  algorithms.  The  method  has  been  applied  to  the  delay  distribu¬ 
tion  analysis  of  the  Capetanakis  Window  Random-Access  algorithm  and  the  Part-and- 
Try  algorithm  both  under  binary  C-NC  feedback.  The  bounds  developed  for  the 
Capetanakis  algorithm,  are  tight  for  all  arrival  rates  within  the  stability  region  of  the 
algorithm.  For  the  Part-and-Try  algorithm,  however,  the  bounds  are  satisfactory  for  rela¬ 
tively  low  arrival  rates.  The  computational  difficulty  in  obtaining  tight  bounds  for  the 
latter  algorithm,  is  due  to  its  complicated  state  space.  The  techniques  can  be  easily 
applied  to  other  Window  Random-Access  algorithms  whose  operating  characteristics  are 


as  described  in  Section  2. 
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CHAPTER  V 


RANDOM  ACCESS  ALGORITHMS  FOR  TWO-CLUSTER 

PACKET  RADIO  SYSTEMS 


V.l  INTRODUCTION 

In  this  chapter  we  consider  a  two-cluster  packet  radio  network.  In  each  cluster  a  sin¬ 
gle  common  channel  is  used  for  packet  transmissions  by  the  users  in  it.  We  assume  that 
the  users  are  not  necessarily  static.  Consequently,  limited-sensing  random  access  algo¬ 
rithms  are  the  most  appropriate  since  they  require  that  a  user  observes  the  feedback  chan¬ 
nel  from  the  time  he  generates  a  packet  to  the  time  this  packet  is  successfully  transmitted. 

In  a  multi-cluster  packet  radio  network  neighboring  clusters  may  overlap.  The  users 
located  in  the  overlapping  regions  are  then  exposed  to  transmissions  and  feedbacks  from 
more  than  one  clusters.  The  latter  phenomenon  can  be  exploited  to  improve  the  delays 
for  the  users  located  in  the  overlapping  regions,  (see  [361  for  a  similar  problem  in  mobile 
telephone  systems).  Consider  for  example  clusters  1  and  2,  and  call  the  users  located  in 
their  overlapping  region,  marginal  users,  let  us  also  call  the  users  in  cluster  i  ,  i  =  1  ,2, 
which  are  not  located  in  the  overlapping  region,  local  users.  The  local  users  in  each  clus¬ 
ter  transmit  their  packets  over  the  common  channel  dedicated  to  that  cluster.  To  resolve 
possible  collisions  they  employ  a  random  access  algorithm,  called  RAAi  (i  stands  for 
cluster  i).  Due  to  the  double  exposure  the  marginal  users  have  a  choice:  they  can  join 
RAAI  or  RAA2  for  the  transmission  of  their  packets.  This  choice  can  be  implemented  in 
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a  static  or  dynamic  way.  The  static  implementation  consists  of  a  priori  assigned  proba¬ 
bilities.  Upon  generation  of  a  new  packet,  each  marginal  user  decides  to  join  RAAi  with 
probability  p,,  where  Pt  +  P2  =  1,  and  remains  there  until  his  packet  is  successfully 
transmitted.  The  policy  described  above  is  simple.  However,  it  requires  that  the  marginal 
users  know  the  a  priori  assigned  probabilities  at  all  times.  The  latter  is  not  considered  a 
serious  disadvantage,  unless  the  a  priori  assignment  changes  with  time  due  to  changes  in 
system  dynamics.  This  is  the  case  of  multi-cluster  systems  with  mobile  users,  where 
users  enter  and  leave  the  system  freely.  In  such  a  dynamically  changing  topology  it  is 
hard  to  update  all  the  marginal  users  regarding  the  optimal  values  of  the  probabilities  p;. 
Consequently  a  dynamic  implementation  should  be  adopted.  Such  an  implementation 
may  in  addition  give  a  delay  advantage  to  the  marginal  users.  The  latter  may  be  desirable 
in  systems  where  the  marginal  users  play  a  special  role,  i.e.  they  transmit  priority  or  con¬ 
trol  messages.  In  this  chapter,  we  consider  two  overlapping  clusters,  each  employing  the 
limited  sensing  RAA  presented  in  chapter  in.  Next,  we  propose  a  dynamic  protocol  via 
v/hich  the  marginal  users  decide  which  cluster  to  join.  The  organization  of  this  chapter  is 
as  follows:  In  section  V.2  we  present  the  system  model.  In  section  V.3  we  describe  the 
protocol  used  by  the  marginal  users.  In  sections  V.4  and  V.5  we  present  the  stability 
analysis.  Section  V.6  contains  our  conclusions. 

V.2  SYSTEM  MODEL 

Consider  a  two-cluster  packet-radio  system.  We  assume  that  in  each  of  the  two 
clusters,  the  synchronous  limited  sensing  random  access  algorithm  presented  in  chapter 
III  is  employed.  Time  is  divided  into  slots  of  length  equal  to  the  packet  duration,  and  the 
starring  instants  of  the  slots  are  identical  in  both  dusters.  By  the  end  of  each  slot  all  the 
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users  in  each  cluster  receive  feedback  information  that  informs  them  about  the  transmis¬ 
sion  activity  in  the  cluster.  This  feedback  is  either  ternary,  (Collision  (C)  versus  Success 
(S)  versus  Empty  (E)),  or  binary,  (Collision  (C)  versus  Non-Collision  (NC)).  In  each 
cluster,  a  local  user  is  required  to  monitor  the  feedback  channel  that  corresponds  to  its 
cluster.  We  assume  negligible  propagation  delays  and  error-free  forward  and  feedback 
channels. 

Each  marginal  user  is  able  to  monitor  both  the  feedback  channels.  At  the  time  when 
a  marginal  user  generates  a  new  packet,  he  starts  monitoring  the  feedback  channels  from 
both  clusters  until  he  decides  to  join  the  operations  of  one  of  the  two  RAAs,  for  the 
transmission  of  his  packet.  Upon  this  decision,  he  maintains  the  continuous  monitoring 
of  the  feedback  channel  that  corresponds  to  the  RAA  he  chose,  until  his  packet  is  suc¬ 
cessfully  transmitted. 

The  mobility  of  the  users  is  assumed  low  so  that,  each  user  remains  wi'hin  the  same 
geographical  region  (local  or  marginal)  from  the  time  he  generates  a  packet  to  the  time 
that  this  packet  is  successfully  transmitted.  Assuming  that  each  of  the  RAAs  operates 
away  from  the  saturation  point  (throughput),  then  with  high  probability  this  time  period 
is  relatively  small. 

It  is  assumed  that  the  local  traffic  generated  in  cluster  i,  i=l,  2,  is  Poisson  distri¬ 
buted  with  intensity  i=l,2,  and  that  the  traffic  generated  by  the  marginal  users  is  also 
Poisson  distributed  with  intensity  X3. 

V.3  THE  ALGORITHM 

We  assume  that  the  two  clusters  employ  the  same  random  access  algorithm.  We 
decided  to  adopt  the  algorithm  presented  in  Chapter  III  since  it  has  simple  operational 
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characteristics  and  is  very  insensitive  to  feedback  channel  errors.  Recall  that  for  the  algo¬ 
rithm  in  Chapter  in,  a  user  with  a  new  packet  waits  for  two  consecutive  noncollision 
(NC)  slots  before  he  transmits  his  packet  for  the  first  time. 

Upon  generation  of  a  new  packet,  a  marginal  user  monitors  both  the  feedback  chan¬ 
nels  until  the  first  time  that  he  is  allowed  to  transmit  his  packet  over  one  of  the  forward 
channels,  (i.e.,  occurrence  of  two  consecutive  NC  slots  in  the  corresponding  feedback 
channel.  If  the  user  is  allowed  to  transmit  his  packet  over  both  the  forward  channels  at 
the  same  time,  he  then  selects  one  of  them  with  probability  0.5).  If  the  received  feedback 
is  NC  then,  the  packet  has  been  successfully  transmitted  and  the  user  stops  observing 
both  the  feedback  channels.  However,  if  the  received  feedback  is  C  then,  the  user  only 
observes  the  feedback  channel  that  corresponds  to  the  forward  channel  over  which  he 
made  his  first  transmission  attempt  and  follows  the  steps  of  the  RAA.  Upon  t^e  success¬ 
ful  transmission  of  his  packet  the  user  stops  monitoring  the  aforementioned  feedback 
channel. 

Upon  generation  of  a  new  packet,  a  local  user  monitors  the  feedback  channel  that 
corresponds  to  its  cluster  until  the  successful  transmission  of  his  packet.  From  the 
description  of  the  first  time  transmission  rule  for  both  the  local  and  marginal  users,  it  is 
clear  that  the  marginal  users  have  an  advantage  over  the  local.  In  particular,  their  waiting 
delay  (time  interval  from  the  packet’s  generation  time  until  the  first  transmission 
attempt)  is  smaller  than  (sometimes  equal  to)  that  of  the  local  users.  Consequently,  their 
total  delays  are  smaller  than  those  for  the  local  users. 
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V.4  ALGORITHMIC  ANALYSIS 

For  notational  convenience  we  will  refer  to  the  local  users  in  cluster  1,  the  local 
users  in  cluster  2,  and  the  marginal  users,  as  subsystem  1,  subsystem  2,  and  subsystem  3, 
respectively.  We  assume  that  the  three  subsystem  traffics  are  generated  by  independent 
processes,  and  that  the  number  of  packet  arrivals  per  slot  in  subsystem  j,  j=  1,2,3,  is  Pois¬ 
son  distributed  with  intensity  Xj.  From  the  throughput  analysis  of  the  adopted  RAA  (see 
Chapter  EH),  we  know  that  the  algorithm  attains  a  maximum  throughput  0.429,  and  that 
the  optimal  window  size  is  A*  =  2.33.  Therefore,  the  following  conditions  are  necessary 
for  stability: 

X,  <  0.429  ,  and  X2  <  0.429  ,  and  X,  +  X2  +  h  <  2  (0.429)  =  0.858  (V.  1 ) 

The  necessary  conditions  in  (V.l)  determine  a  {Xj,  j=l ,2, 3 }  hyperplane,  which  con¬ 
tains  the  {Xj,  j=l,2,3J  region  that  determines  the  system  throughput.  Bounds  on  the 
(Xj,  j=l,2,3}  space  which  provides  the  system  throughput,  will  be  obtained  via  the  sta¬ 
bility  analysis  in  Section  V.5  below. 

V.5  SYSTEM  STABILITY 

Assume  that  the  algorithms  in  both  clusters  start  operating  at  time  zero.  Consider 
the  sequence  in  time  of  the  collision  resolution  intervals  (CRIs),  induced  by  the  two 
RAAs.  We  define  the  sequence  {T,,},^  as  follows:  (1)  Tn  corresponds  to  the  starting 
slot  of  some  CRI,  (  note  that  at  Tn,  two  CRIs  may  simultaneously  begin;  one  for  each  of 
the  two  RAAs).  (2)  To=2,  and  at  T0  two  CRIs  begin;  one  for  each  RAA. 

Let  {T^ *  Jnso  be  the  subsequence  of  the  sequence  {Tn)n>o,  which  consists  of  those 
time  instants  when  two  CRIs  begin  simultaneously.  Clearly,  T{>5)=To=2.  Let 
j= 1,2,3,  denote  the  total  length  of  the  unresolved  arrival  intervals  in  subsystem  j,  at  the 
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time  instant  T^s).  D^s  is  then  called  "the  lag  of  subsystem  j  at  time  T^s)."  From  the 
operations  of  the  algorithm  we  conclude  that:  (1)  D^s>2  and  that  the  possible  values  of 
D^s,  j  =1,2  are  countably  many.  (2)  D$s=2,  j=l,2,3.  (3)  At  time  T^s),  the  RAA  in  clus¬ 
ter  k,  k=l,2,  examines  two  arrival  intervals:  one  from  subsystem  k  which  has  length 
min(D^,  A)  and  contains  arrivals  generated  by  a  Poisson  process  with  intensity  Xk 
arrivals  per  slot,  and  one  from  subsystem  3  which  has  length  min  (D^s>  A)  and  contains 
arrivals  generated  by  a  Poisson  process  with  intensity  0.5  X3  arrivals  per  slot.  (4)  The  tri¬ 
ple  (D^s,  j=l,2,3)  describes  the  state  of  the  system  at  time  and  the  sequence 
{Sn)n>o^  {D^,j=l,2.3}n>o  is  a  three-dimensional  irreducible  and  aperiodic  Markov 
Chain. 

The  stability  of  the  system  is  determined  by  the  ergodicity  of  the  three-dimensional 
Markov  Chain  {S,,},^.  In  addition,  at  time  T(„s\  the  backlogs  of  each  of  the  three  sub¬ 
systems  are  represented  by  the  three  lags,  D^s,  j=l,2,3.  Given  D^s  =d$\  j=l ,2,3,  we 


define  the  expected  system  backlog  at  time  T^s)  to  be  the  ]£  Xjd$\  where  V(dn)  = 


V((d^}iSjS3)  _  jr  Xj  d^  is  a  Lyapunov  function  of  the  three  subsystem  lags  (d^ }  ^<3. 
j=i 


Let  C  denote  the  state  space  of  the  Markov  Chain  (Sn  }„><),  and  define  the  operator 
AV(d)  =  AV^d^}]^),  called  a  generalized  drift,  as  follows: 


AV(d)  =  AV({d<J))1<J<3)  =  E{V({D^1,s)1<J<3)-V({DW}1<J<3)lD^=d^,j=l,2,3} 

-E(£xjD^I.sID<i>s=dW,j=l,2,3}-ixjd^;d={d^}I<j<3  e  C  (V.2) 

j=i  i=i 

Since  we  are  interested  in  determining  the  system’s  throughput  we  assume  that  the  sys¬ 
tem  operates  at  the  saturation  point  therefore,  all  the  lags  are  assumed  to  be  sufficiently 
long,  so  that  in  the  time  interval  [Tns),  T^j],  the  examined  interval  is  always  of  length  A. 
Consequently  the  number  of  packets  that  are  successfully  transmitted  during  the  first  CRI 
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that  belongs  to  RAAi,  i  =  1,2  in  the  above  interval,  is  Poisson  distributed  with  parame¬ 
ter  (  Xi  +  0.5  A-3  )  A.  Furthermore,  the  number  of  packets  that  are  successfully  transmit¬ 
ted  during  every  subsequent  CRI  that  belong  to  RAAi,  i  =  1,2  in  the  above  interval,  is 
also  Poisson  distributed  with  parameter  (  Xi  +  X3  )  A.  Under  the  above  conditions  and 
given  the  subsystem  Poisson  rates  {Xj }  i<j<3 ,  let  us  define: 

Ex{/  I  u}:  Given  that  the  length  of  the  interval  to  be  examined 

equals  to  u  and,  that  the  Poisson  input  rate  equals  X, 
the  expected  number  of  slots  needed  for  the  suc¬ 
cessful  transmission  of  all  the  packet  arrivals  within 
the  examined  interval. 

i=l,2  ;  ),A):  The  number  of  CRIs  in  [Tf,s),  that  belong  to 

RAAi,  i  =  1,2. 

AV(dA,  {A.j }):  The  generalized  drift  in  (V.2).  We  next  define  a 

generalized  drift  for  the  RAA  in  the  i-th  cluster: 

i=l,2;  A(i)V(dA,  {Xj})^  Ex^.2  (/I A}-A  + 

+  [E{Nl((Xj},A)}-l]  [e^+x,  { / 1 A } —  aJ  (V.3) 

Then,  from  (V.2)  we  derive  the  following  expression: 

—  \  \ 

AV(dA,  {Xj})  =  (X.1+-y-)A<1>V(dA,  {?iJ})  +  (X2+-y-)A^V(dA,{XJ)H 

- ^~^E{N,({A.j),A)}  +E(N2([Xj  ),  A) ) -2]  (V.4) 

where, 


ExI+x,i2{/IA)+[E{N1({XJ),A))-l]EX|+Xj{/IA}= 
=  EXi+x,i2{/IA)+[E{N2({Xj),A))-1]EX2+X3{/IA) 


(V.5) 
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We  now  state  a  Lemma  whose  proof  is  similar  to  the  one  presented  in  [37], 

Lemma  1 

(i)  Let  there  exist  some  £>0,  such  that  the  two  conditions  below  are  satisfied: 

A(I)V(dA,{Aj])  <-e 

A(2)V(dAt  {Xj))<-e  (V.6) 

Then,  the  Markov  chain  {Sn  }n>o  is  ergodic  at  the  Poisson  rates  { A-j )  i<j<3  . 

(ii)  Let  at  least  one  of  the  two  generalized  drifts,  AwV(dA,  (Aj}),i=l,2,  be  nonnega¬ 
tive.  Then,  the  Markov  Chain  {Sn}n>o  is  nonergodic  at  the  Poisson  rates 
{^jh<j<3  •  O 

The  conditions  in  (V.6)  define  lower  bounds  on  the  |  Aj )  regions  for  which  the  Mar¬ 
kov  chain  {Sn  )n>o  is  ergodic. 

In  Appendix  D,  we  describe  the  methodology  used  in  the  evaluation  of  the  max¬ 
imum  (Aj)  values  for  which  the  chain  is  ergodic.  In  Figure  V.l,  we  plot  the 

boundaries  of  the  (A(,  A2)  stable  regions,  parametrized  by  various  A3  values.  Those 
boundaries  are  symmetric  around  the  45°  straight  line.  In  Figure  V.2,  we  plot  A3  against 
A  =  At  =  A2,  (symmetrically-loaded  system). 

From  Figures  V.l  and  V.2,  we  observe  that  the  maximum  value  of  the  sum 
A1+A2+A3  is  always  strictly  less  than  0.858  (i.e.,  strictly  less  than  twice  the  throughput  of 
each  local  RAA).  This  is  so  because  under  the  conditions  described  above  the  RAAi, 
i  =  1,2,  operates  with  an  input  rate  equal  to  (  A^  +  A3  )  A  during  most  of  the  time.  Conse¬ 
quently,  if  the  rate  of  the  local  traffic  is  close  to  the  algorithmic  throughput  for  one  of  the 
RAAs  in  the  system,  the  introduction  of  marginal  users  results  in  instability. 
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V.6  CONCLUSIONS 

In  this  chapter,  we  performed  stability  analysis  for  a  two-cluster  packet  radio  net¬ 
work.  Each  cluster  contains  local  and  marginal  users.  In  each  cluster  a  single  common 
channel  is  available  for  packet  transmission.  Packets  that  are  generated  by  the  local  users 
in  a  cluster,  can  be  only  transmitted  over  the  forward  channel  that  corresponds  to  that 
cluster.  In  contrast,  packets  that  are  generated  by  the  marginal  users  can  be  transmitted 
over  either  one  of  the  two  available  forward  channels.  We  propose  a  protocol  according 
to  which  the  marginal  users  are  able  to  dynamically  select  the  forward  channel  over 
which  they  will  eventually  transmit  their  packet.  According  to  the  protocol  rules,  each 
marginal  user  with  a  packet  for  transmission  is  required  to  observe  both  the  feedback 
channels  from  the  packet’s  generation  time,  in  order  to  decide  over  which  channel  he  will 
transmit  his  packet.  No  a  priori  knowledge  of  the  input  traffic  rates,  or  the  subsystem’s 
states  is  required. 


l 

1 

I 


82 


REFERENCES 


[1]  D.  Bertsekas,  R.  Gallager,  Data  Networks.  Prentice  Hall,  1987. 

[2]  P.  E.  Jackson  and  C.  D.  Stubbs,  "  A  Study  of  Multiaccess  Computer  Communica¬ 
tion,"  in  1969  Spring  Joint  Comp.  Conf.  AFIPS  Conf.  Proc.,  Vol.  34,  Monivale  N.J, 
AFIPS  Press,  pp.  491-504. 

[3]  N.  Abramson,  "  Packet  Switching  With  Satellites,"  in  1973  Nat.  Comp.  Conf.,  AFIPS 
Conf.  Proc.,  Vol.  42,  New  York,  AFIPS  Press,  1973,  pp.  695-702. 

[4]  N.  Abramson,  "  The  ALOHA  System  -  Another  Alternative  for  Computer  Communi¬ 
cations,"  1970  Fall  Joint  Comp.  Conf.,  AFIPS  Conf.  Proc.,  Vol.  37,  pp.  281-285. 

[5]  G.  Fayolle  and  E.  Gelenbe,"  Stability  and  Optimal  Control  of  the  Packet  Switching 
Broadcast  Channel,"  Journal  of  the  Association  for  Computing  Machinery,  Vol.  24,  No. 
3,  July  1977,  pp.  375-386. 

[6]  B.  S.  Tsybakov  and  V.  A.  Mikhailov,  "  Ergodicity  of  a  Slotted  ALOHA  System," 
Problemy  Peredachi  Informatsii,  Oct.-Dec.  1979,  Vol.  15,  pp. 73-87. 

[7]  J.  I.  Capetanakis,  "Generalized  TDMA:  The  Multiple  Accessing  Tree  Protocol,"  IEEE 
Trans,  on  Communications,  Oct.  1979,  COM-27,  pp.  1476- 1484. 

[8]  B.  S.  Tsybakov  and  V.  A.  Mikhailov,  "  Free  Synchronous  Packet  Access  in  a  Broad¬ 
cast  Channel  With  Feedback,"  Problemy  Peredachi  Informatsii,  Oct.-Dec.  1978,  Vol.  14, 
No.  4,  pp.  34-59. 

[9]  R.  G.  Gallager,  "  A  Perspective  of  Multiaccess  Channels,"  IEEE  Trans,  on  Inform. 
Theory,  Special  Issue  on  Random  Access  Communications,  March  1985,  Vol.  IT-31, 


83 


[10]  B.  S.  Tsybakov  and  V.  A.  Mikhailov,  "  Random  Multiple  Access  of  Packets,  Part- 
and-Try  Algorithm,"  Problemy  Peredachi  Informatsii,  Oct.-Dec.  1980,  Vol.  16.  No.  3,  pp. 
65-79. 

[11]  J.  L.  Massey,  Collision  Resolution  Algorithms  and  Random  Access  Communica¬ 
tions,  in  Multi-User  Communication  Systems,  Ed.  G.  Longo,  CISM  Coupes  and  Lec¬ 
tures,  No.  265,  pp.  73-137,  Springer-Verlag,  New  York,  1981. 

[12]  B.  S.  Tsybakov  and  N.  D.  Vvedenskaya,  "  Random  Multiple  Access  Stack  Algo¬ 
rithm,"  Problemy  Peredachi  Informatsii,  July-Sept.  1980,  Vol.  16,  No.  3,  pp.  80-94. 

[13]  N.  D.  Vvedenskaya  and  B.  S.  Tsybakov,  "  Random  Multiple  Access  of  Packets  to  a 
Channel  With  Errors,"  Problemy  Peredachi  Informatsii.  1983,  Vol.  19,  No.  2,  pp.  52-68. 

[14]  P.  Mathys,  Analysis  of  Random  Access  Protocols,  Ph.D  Thesis,  Swiss  Federal  Insti¬ 
tute  of  Technology,  Zurich,  1984. 

[15]  L.  Merakos  and  C.  Bisdikian,  "  Delay  Analysis  of  the  n-ary  Stack  Algorithm  for  a 
Random  Access  Broadcast  Channel,"  Proceedings  of  the  22nd  Allerton  Conf.  on  Comm. 
Comp.,  Univ.  of  Illinois  at  Urbana-Champaign,  Illinois,  October  1984,  pp.  385-394. 

[16]  L.  Georgiadis,  Limited  Sensing  Random  Access  Algorithms  ana  Unified  Methods 
for  their  Analysis,  Ph.D  Dissertation,  Univ.  of  Connecticut,  Storrs.  1986. 

[17]  L.  Georgiadis  and  P.  Papantoni  -  Kazakos,  "  A  0.487  Throughput  Limited  Sensing 
Algorithm,"  IEEE  Trans,  on  Inform.  Theory,  Vol.  IT-33,  No.  2,  March  1987. 

[18]  P.  A.  Humblet,  "  On  the  Throughput  of  Channel  Access  Algorithms  With  Limited 
Sensing,"  IEEE  Trans,  on  Communications,  Vol.  COM-34,  No.  4.  April  1986. 

[19]  B.  S.  Tsybakov,  "  Survey  of  USSR  Contributions  in  Random  Multiple  Access  Com¬ 
munications,”  IEEE  Trans,  on  Info.  Theory,  Vol.  IT-31,  No.  2,  March  1985. 

[20]  L.  Georgiadis,  L.  Merakos,  and  P.  Papantoni  -  Kazakos.  "  A  Method  for  Delay 


84 


IEEE  Journal  on  Selected  Areas  in  Communications,  July  1987. 

[21]  J.  Kurose,  M.  Schwartz,  and  Y.  Yemini,  "  Multiple  Access  Protocols  and  Time  Con¬ 
strained  Communication,"  ACM  Computing  Surveys,  Vol.  16,  No.  1,  March  1984. 

[22]  J.  Kurose,  M.  Schwartz,  and  Y.  Yemini,  "  Controlling  Time  Window  Protocols  for 
Time  Constrained  Communication  in  a  Multiple  Access  Environment,"  in  Proceedings  of 
the  8th  International  Data  Communications  Symposium,  Oct.  1983,  pp.  75-84. 

[23]  B.  S.  Tsybakov,  N.  B.  Likhanov,  "  Upper  Bound  for  the  Delay  in  a  Random  Access 
System  With  a  Splitting  Algorithm,"  Problemy  Peredachi  Informatsii,  Vol.  18,  No.  4,  pp. 
76-84,  Oct.-Dec.  1982. 

[24]  J.  C.  Huang  and  T.  Berger,  "  Delay  Analysis  of  Interval  Searching  Contention  Reso¬ 
lution  Algorithms,"  IEEE  Trans,  on  Inform.  Theory,  Vol.  IT-31,  No.  2,  March  1985. 

[25]  G.  Polyzos,  M.  Molle,  and  A.  Venetsanopoulos,  "  Performance  Analysis  of  Finite 
Non-Homogeneous  Population  Tree  Conflict  Resolution  Algorithms  Using  Constant  Size 
Window  Access,"  IEEE  Trans,  on  Comm.,  Vol.  COM-35,  No.  1,  pp.  1124-1138, 
November  1987. 

[26]  J.  I.  Capetanakis,  "Tree  Algorithms  for  Packet  Broadcast  Channel,"  IEEE  Trans. 
Inform.  Theory,  Vol.  IT-25,  pp.  505-515,  Sept.  1979. 

[27]  R.  G.  Gallager,  "Conflict  Resolution  in  Random  Access  Broadcast  Networks",  in 
Proc.  AFOSR  Workshop  Communication  Theory  and  Applications,  Provincetown,  MA, 
Sept.  1978,  pp.  74-76. 

[28]  M.  Paterakis,  L.  Georgiadis,  and  P.  Papantoni-Kazakos,  "  On  the  Relation  Between 
the  Finite  and  the  Infinite  Population  Models  for  a  Class  of  RAAs",  IEEE  Trans,  on 
Comm.,  Vol.  COM-35,  pp.  1239-1240,  Nov.  1987. 

[29]  P.  Studer,  and  H.  Pletscher,  "Q-ary  Part-and-Try  Algorithm  for  Packet  Conflict 
Resolution",  Inst.  Telecommun.  Swiss  Fed.  Inst.  Technol.,  Zurich,  Switzerland,  Oct.  1 1, 


85 


1984. 

[30]  A.  G.  Pakes,  "Some  Conditions  for  Ergodicity  and  Recurrence  of  Markov  Chains", 
Oper.  Res.,  Vol.  21,  pp.  617-622,  1973. 

[31]  M.  Kaplan,  "A  Sufficient  Condition  for  Nonergodicity  of  a  Markov  Chain",  IEEE 
Trans.  Info.  Theory,  Vol.  IT-25,  pp.  470-471,  1979. 

[32]  T.  Gonsalves,  "Packet-Voice  Communication  on  an  Ethernet  Local  Computer  Net¬ 
work:  An  Experimental  Study,"  In  Proc.  of  the  SIGCOM  Communications  Architectures 
and  Protocols  Symposium,  (Austin,  TX,  March  1983).  ACM,  New  York.  pp.  178-185. 

[33]  G.  Nutt,  and  D.  Bayer,  "Performance  of  CSMA/CD  Networks  Under  Combined 
Voice  and  Data  Loads,"  IEEE  Trans,  on  Comm.,  Vol.  COM- 30,  Jan.  1982,  pp.  6-11. 

[34]  E.  Seneta,  Non-negative  Matrices.  John  Wiley&Sons,  New  York,  1973. 

[35]  Ayal  Bar-David,  and  Moshe-Sidi,  "Collision  Resolution  Algorithms  in  Multistation 
Packet  Radio  Networks,"  submitted  for  publication. 

[36]  B.  Eklundh,  "Channel  Utilization  and  Blocking  Probability  in  a  Cellular  Mobile 
Telephone  System  with  Directed  Retry,"  IEEE  Trans.  Comm.,  Vol.  COM-34,  pp.  329- 
337,  April  1986. 

[37]  W.  Szpankowski,  "Stability  Conditions  for  Multidimensional  Queueing  Systems 
with  Applications,"  J.  of  Operations  Research,  to  appear. 

[38]  W.  Szpankowski,  "Some  Theorems  on  Instability  with  Applications  to  Multi-Access 
Protocols,"  J.  of  Operations  Research,  to  appear. 

[39]  R.  L.  Tweedie,  "Criteria  for  Classifying  General  Markov  Chains,"  Adv.  Appl.  Prob., 
Vol.  8,  pp.  737-771,  1976. 


I 


$6 


APPENDIX  A 

A.1  Recursions 

Towards  the  computation  of  E{Qd }  and  P(l\u,d),  let  us  denote  by  E{\kl,k2,B} 
the  expected  value,  given  that  £;  packets  have  counter  value  equal  to  i,  i=l,2,  and  that  the 
maximum  length  of  the  collision  resolution  interval  is  B.  Let  P  {-\k\,k2,B)  denotes  con¬ 
ditional  probability,  given  i=0, 1,  and  B  as  above.  Then,  for  t \ud  denoting  either  one 
of  the  quantities  luj,  zUid,  and  nu  d,  we  clearly  have: 

E(Ud)=Z  E{^d\k,0J~d}e~u 

k- 0  K  ■ 

P  (/ 1  u,d)  =  £  P  (/ 1  k,  0  J-d)e~^ 

k= 0  k  ' 

It  is  also  easy  to  prove  that, 

E(Vu,d}  =  E{nu  d} 

Note  that  sup  E{C,Utd  l/fc,  0,7 -d}  <  Therefore,  the  computation  of  a  finite  number 

of  terms  is  sufficient  for  developing  tight  bounds  on  the  quantities  of  interest.  Also,  note 
that  given  k\,  k2,  B,  the  operation  of  the  algorithm  is  independent  of  u  and  d.  Let  us 
denote: 


Ek\,ki,B  ~  E{ludlk  i  ,k2,B},  E{nudlk\,k2,B},  Z^^^^g  —  E{zudlk  \,k2,B} 

From  the  operation  of  the  algorithm  we  derive  the  following  recursions: 
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Lk„k2.B  = 


1  +  *^.0,8-1  .  0<ki<l 


1  +Li,k1+ki-i.B-l  »  W-P- 


2"\  *,>1, 


*  j  >/>0 


tyt„*2,8 


*42,o.b-i,  *i=0 

***2,0,B-1  +  1.  *1  =  1 

r  > 


**«.*, +*2-1.8 -1« 


w.p. 


,  Jtj  >1,  *!>i>0 


1+***2.0.B-1  +Z*2,0,B-1.  *i  =  l 

**t2, 0,8-1  +  Z*2. 0,8-1-  *1=0 


**i.*,+*2-i.B- 1  +  Zi,*t+*2-i,B-l.  W.p. 


r  *\ 
* 


2  *',*!>!, *!>2>0 


/J(/I*1,*2,fl)=^ 


P(/-ll*2, 0,5-1),  k  i=0  or  1 
/*  (/ — 1  Ij,*i+*2-/,£-1),  w.p. 


y~k 


2  1  ,  *  |  >  1 ,  *  1  >/  >0 


The  initial  conditions  in  the  above  recursions  are  the  following: 


*-0,0 ,B  =  *>1 ,0,B  =  *'2,0,1  =*'0.1.1  =  **-l 

**2,0,2  =  *-0,1,8  =  2,  B> 2 
**0,0,8  =**0,1,1  =**2,0,1  =0’  **-l 
**i.o.b  =  1,  5>1,  **2,o,2  =  0-5,  5>2 
**0.1.8  =  1,5>2 
Z 1.1.1  =  Z  1,0,8  =  1.  **^1 
Z*,,*2.i  =0,  *t>2,  *2^0,  Zij.ii  =  3,  B> 2 

Zk,.k2.2  =  2* , 2-*1 ,  k\>2 

P',l\kl,kz,B)  =  <:,  1>B 
P( 2  I  1,0,5)  =  P(/ 1  1,0,5),  5>1,  /> 2 
5(11*,  0,5)  =  0,  5>2,  *>2 


I 
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/>(lll,0,fl)  =  />(llifc,  0, 1)  =  P(2\k,  0,2)=  1,  B>1,  k> 2 
P(2\k\,Jc2,2)=  1,  ki>2 

A.2  Computations  in  the  Absence  of  Bounds  on  Delays 

Let  T  =°o,  and  let  us  then  define: 


Lruk_n:  The  expected  number  of  slots  needed  by  the  algorithm,  for  the  successful 
transmission  of  k  packets,  given  that  n  of  the  packets  have  counter  values  equal  to  one, 
and  k-n  of  the  packets  have  counter  values  equal  to  two. 

From  the  operation  of  the  algorithm,  we  then  obtain: 

L/i.i-n  —  1  +  £ I" ^  2  "L i.jk— i  (A.l.a) 

i=0L  J 

Lo,o  =  Liio  =  L  Lo.,  =  1+L,  q;  1^1  (A.l.b) 

It  can  be  shown  by  induction  that  has  the  following  form: 

Ln.t-n-A^Lt.o+^Lt.j.o+A?);  *22,  2 <n<k  (A.2) 

;  where  A%\  i  =  1,2,3  are  independent  of  k  and  can  be  computed  recursively  as  fol¬ 
lows: 


A^m-2-r1  ^jl+sV’flj,  n>3  (A. 3. a) 

A^l-2-"]-1  2-”| «  +  IlA<2)|^j,  n>3  (A.3.b) 

A^3>=[l-2-n|l+2-"(n+l)+2-',',£  AP>[^|,  n>3  (A.3.C) 

A^=  1/3,  Ai2)  =2/3,  Ap  =7/3 


(A.3.d) 
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For  n=k,  expression  (A.2)  gives: 


E{lu.<i}=iUe-h‘^-  (A.7) 

t=o  K- 

In  the  same  way,  we  can  develop  recursive  formulas  and  bounds  for  the  quantities 
associated  with  the  cumulative  delays  of  the  packets  during  a  CRI.  These  quantities  can 
then  be  used  for  the  computation  of  the  delays  in  the  unconstrained  case,  by  following 
the  methodology  presented  in  [20]. 
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APPENDIX  B 


B.l.  Stability  Analysis 

We  present  the  stability  analysis  in  the  case  feedback  errors  may  occur.  By  setting 
e  =  8  =  0  in  our  results,  we  get  the  corresponding  quantities  in  the  error  free  case. 

The  stability  region  of  the  algorithm  is  provided  by  inequality  (131.1),  where  the 
expected  value,  is  given  by  expression  (III.2).  We  start  with  the  computation  of  the 
expected  values  L*  ;  k  >  0  in  (III. 2a). 

Computation  of  Lt 

We  define: 

Gn.k-n-  The  expected  number  of  slots  needed  by  the  algorithm,  for  the  successful 
transmission  of  k  packets,  given  that  n  of  those  packets  have  counter  values  equal  to 
one,  and  k-n  of  the  packets  have  counter  values  equal  to  two. 

Notice  that  Lk  =  G^o-  for  k> 2,  while  L *  *  G^o-  for  k  <2.  We  first  show  how  to 
compute  Lq  and  L  i . 

Computation  of  Ln\ 

From  the  operation  of  the  algorithm  we  have 


where 


J  1  ;  w.p.  (1-e) 

1  1  +  G o.o  ;  w.p.  e 


(B.  1.1) 
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Go.o  - 


\  +Lq  4,  w.p.  (1— e) 
i  +  Gq.o  ;  wP-  e 


From  (B.1.1)  and  (B.1.2)  we  find  that 


L  o  = 


1 

(1-e)2 


Computarion  of  L 


We  find  that  L ;  satisfies  the  following 


1  ;  w.p.  (1-5) 

L  i  =■  1+G  u  ;  w.p.  0.55 
1+Go.i  ;  w.p.  0.55 

where  G  i  o  and  Go. i  satisfy  the  following 


1+L0  ;  w.p.  (1-5) 
G  i,o  1 +G  i,o  ;  w.p.  0.55 

1+G  o,i  ;  w.p.  0.55 

r  i+£>i  i  w.p.  (i— e) 

G°a  =  ]  1+G 0,1  ;  w.p.  e 


From  (B.1.4),  (B.1.5),  and  (B.1.6)  we  find  that 

-l 


L,= 


1- 


2-5 


1+ 


2(l-e)  (2-5) 


1+ 


(1-5)  ,  5 

(1-e)2  2(1  e) 


Computation  of  Lv,  for  k>2. 

From  the  operation  of  the  algorithm  we  obtain: 


where 


Gn.k-n  =  1  +  Z 

i  =0 


2~n  G,_k~i 


n> 2  ,  k>n 


Go .k  - 


1 

(1-e) 


+  Lk 


(B.1.2) 

(B.  1.3) 

(B.1.4) 

(B.1.5) 

(B.1.6) 

(B.1.7) 

(B.  1.8) 


(B.1.9) 
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G  it  = 


l+L*(l-6)  +  0.58  - +  Z.£+ 1 

(*  e) 


(B.I.10) 


1-0.58 

It  can  be  shown  by  induction  that  G^k-n  has  the  following  form 


G^-n  =  ^1)Gfc0+^2)GJt_1.0+A(l3)  ;  2<n<k  (B.1.11) 

;  where  A i=l ,2,3  are  independent  of  k  and  can  be  computed  recursively  as  follows: 


4m  2+5  n)  _  4(1-8)  A(3)  _  14-35- 12e  +  4e5 

*  '  6-35  ’  2  6-35  *  2  3(1— e)  (2-5) 


^>  =  [I-2-T>  2-  1^  + 


^>  =  n-rv  + 


(B.  1.12) 


(B.  1.13) 


(B.1.14) 


A»,  =  [1-2'"r1{1+  71^6)  +  (5)-+  (2^)  +2"|(/]  ^t’^3  (B.  1.15) 


For  n=k,  expression  (B.1.11)  gives: 


AP  Ai» 


(B.  1.16) 


Expression  (B.  1.16)  together  with  the  recursions  in  (B.1.13)-(B.1.15)  provide  a  mean  for 
the  computation  of  L *,  k>2. 


Development  of  an  upper  bound  on  Li 


It  can  be  seen  by  induction  that: 


4'>S}+2-8 

(B. 1.17) 

a[2) 

-^—=l-,k>2 

1-Alt1  > 

(B.  1.18) 

p  (2-2E+8)  k+  1 

*  (2— 5)(  1— e)  (1-e) 

(B.  1.19) 

From  (B .  1 . 1 6)-(B.  1.19)  we  find: 
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(1 — e)(4— 55)  (l-e)(4-55) 


(B.  1.20) 


from  which  we  can  show  that 


3(2-2e+S)  2  3(6-'2e-8) 

*  “  2(l-e)(4-58)  2(  1  — e)<4— 5S) 


L  i  - 


6(2— e) 


(l-e)(4-58) 


^Lut  =  aJc2^k+y  (B.  1.21) 


Due  to  the  upper  bound  on  Lk,  we  conclude  then  that  the  following  condition  is  sufficient 
for  stability: 


where 


30 


I  Lkp{k  IA)+  X  L\  p(k  I  A)  <  A 

k  = 0  *  = 31 


(B  1.22) 


p(k\A)  =  e  ^ 


(XA)* 

it! 


After  some  manipulations,  we  conclude  that  (B.1.22)  is  equivalent  to: 


30  f  30  1 

f(XA)  =  X  4p(/tlA)  +  a^(XA)2-f  AA-  X  k2p(klA)\  + 

k-0  [  t=0  J 


30  "j 

30 

+  p 

XA-  X  k  pik\ A) 

+  Y 

1-  Xp&  1  A) 

*=o  J 

l  J 

Let  us  now  define: 


(B.1.23) 


x^XA  (B.1.24) 

Then,  from  (B.  1.23)-(B.1.24)  we  conclude  that,  for  the  stability  of  the  algorithm,  it  is  suf¬ 
ficient  that  the  input  rate  X  satisfies  the  following  inequality. 

X<sup  (B.  1 .25) 

x>o  fix) 

The  following  condition  specifies  a  region  of  X  values  for  which  the  algorithm  is 
unstable. 


30 


where  g(.x)  =  X  Lk  e 


~x  xk 


Ik ! 


*=o 


X  >  sup  — — 
*>o  gU) 


(B.  1 .26) 


I 
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l 


The  maximization  of  expessions  in  (B.1.25)-(B.1.26)  has  been  done  numerically,  and 
provides  the  throughput,  as  well  as  the  optimal  window  size  A* .  In  all  the  caseo  the  order 
of  the  difference  between  the  two  suprema  in  (B.1.25)-(B.1.26)  is  less  than  10-3.  The 
optimal  window  size  is  found  asx*  (X*)-1,  where  x  *  is  the  value  that  attains  the  suprema 
in  (B.1.25). 


Computation  of  Lu.  for  the  Cacetanakis’s  Dynamic  Algorithm 


Here,  the  quantitites  Lk,  k>2  can  be  computed  recursively  as  follows: 


Lk 


=  [  1-2  2-*rl  j  1+2  L'o  2~k  +  2  *£  L] 


r*i 


,  k>2 


(B.  1 .27) 


where 


L0 


_1 _ 

d-2e) 


and  L\  - 


1-2e  +  5 
(l-2e)(l-S) 


Moreover  the  following  upper  bound  on  Lk,k> 4,  has  been  found. 


Lk  — 


3(  1— e)  |  2(5— e) 

(l-2e)  (l-2e)(l-6) 


(B.1.28) 


(B.  1 .29) 


Using  the  bounds  in  (B.1.29),  we  computed  upper  and  lower  bounds  on  the  throughput 
for  the  Capetanakis’s  dynamic  algorithm.  The  computed  upper  and  lower  bounds  were 
identical  to  each  other  up  to  the  fourth  decimal  point,  and  are  included  in  Table  III. 3. 


B.2.  Delay  Analysis 

The  following  definitions  will  be  used  in  the  sequel. 

/  :  Length  of  a  collision  resolution  interval  (CRI). 
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x  :  Window  size. 

E(Xlx)  :  Expected  value  of  the  random  variable  X,  given  that  the  window  size  is 

X. 

p  (x/x)  :  The  probability  that  the  CRI  has  length  x  given  that  the  window  size  is  x. 

gd  :  Number  of  slots  needed  to  reach  a  CRP  with  lag  I  given  that  the  current 

lag  is  equal  to  d,  de  F. 

wd  :  Cumulative  delay  experienced  by  all  the  packets  that  were  successfully 

transmitted  during  gd  slots. 

N  :  Number  of  packets  transmitted  in  the  CRI,  that  starts  at  time  t. 

z  :  Cumulative  delay  of  the  N  packets,  after  the  CRP  t. 

\\r  :  Cumulative  delay  of  the  N  packets,  until  the  instant  1 2. 

Let  us  also  define. 


Note  tha'  by  definition. 


Gd  -  E  (gd)  Wd  =  E{wd) 


Qx 

Wl=E('£Dn)  =  W 

n= 1 

Also,  the  mean  session  length  G  =  E  (ri+1-T,) ;  i  >  1,  equals  to  G  ] .  If  G 1  <  °®,  then  by 
Wald’s  identity  we  have  that. 


Q=XG, 

Therefore,  the  determination  of  Wj  and  G]  will  permit  the  computation  of  the  mean 
packet  delay. 

The  operation  of  the  algorithm  yields  the  following  relations  for  the  d  e  F. 


1  <  d  <  A  gd  =  l  if  /  =  1 


(B.2.1) 


I 
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I 


I 


i 


I  <d  £A  gd  =  l  +  gi  if  /  >  1  (B.2.Ia) 

and  that 


A  <d  gd  =  l  +gd-A  +  i  (B.2.2) 

Taking  expectations  in  (B1)-(B2)  yields: 


Gd  =  E  (l  \d)  +  X  GsPisId)  if  l<d<A,  de¥ 

s*\.sef 

Gd  =  E(l\A)  +  ZGd-*  +  sp(s/ A)  if  d>A,  deF 
se  F 


(B.2.3) 

(B.2.4) 


Equations  (B.2.3)  and  (B.2.4)  comprise  a  denumerable  system  of  linear  equations.  Of 
interest  to  us  is  the  element  G  \  of  a  particular  solution  of  this  system.  We  now  proceed 


in  the  development  of  an  initial  upper  bound  on  the  solution  of  the  system  in  (B.2.3)- 
(B.2.4).  Following  the  methodology  in  [20],  such  a  bound  will  be  the  sequence 
Gd.u  =  Y ud+C,u<  if  Yu.  Cu.  can  be  determined  so  that  the  following  inequalities  are  satis¬ 
fied 


G2.„  >£(/!<*)+  I  C?.up(s/d)  =  Ci,u  \  <d  <A,  de  F 

j*l..re  F 

G°du  >£(/  IA)+  XG°d-&  +  ,p(s/A)  =  G]d'U  if  d>A,  de  F 
J6  F 


(B.2.5) 

(B.2.6) 


Substituting  Gd  in  the  right  hand  side  of  inequalities  (B.2.5)  and  (B.2.6),  it  can  be  seen 
that  if  de  F, 


Gxd,u=Glu  +  iu(E(Ud)-d-(\+'kd)e-u)~^(,\+Xd)e-Xji  if  1  <d<A  (B.2.7) 

Gld'U=G°d.u+E(l/A)-yu(A-E(l/A ))  i fd>A  (B.2.8) 


From  (B.2.8)  we  conclude  that  if  E  (//A)  <  A,  the  condition  for  stability  of  the  system, 
inequalities  (B.2.6)  are  satisfied  if, 


£  (/ 1  A) 

Y“~A-£(//A) 


(B.2.9) 


With  this  value  of  yu,  it  can  be  seen  that  inequalities  (B.2.5)  are  satisfied  if 
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where. 
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E(z  \d)  +  E(y\d)  +  ^(EO1  \d)  -dl-(  1  +  \d)e~“)  -  vu(d  -E{1 1 d)  -  (1  +  ld)e 

m- - - 

Substituting  Wd%  in  the  right  hand  side  of  (B.2.15MB.2 .16),  we  obtain  Wd_ .  We 


W\u  =  E(z  1 1)  +  £(yl  1)  +  1 1H1  +  \)e^)  +  vH(E(l  1 1)  -  (1  +  X)e'x)  +  ^u(l  -  (1  +  \)e~x ) 


W'u  =  <«-(^-^)(i-(i  +  X)e^) 


(B.2.23) 


G\'U  =  E{1  1 1)  +  yu(E(l  1 1)  -  (1  +  X)e~k)  +  £„(!  -  (1  +  X)e~x)  (B.2.24) 


G}./  =  G}.„-(Cll-C/)(l-(l+X)e->') 

From  the  regenerative  theorem  [20],  we  have  that 


(B.2.25) 


.  W  ,  W\u 

d'  =  — -r-  <  D  <  — =  D“ 

XG  „  XG  / 


(B.2.26) 


In  this  Appendix  we  also  show  that  the  conditional  expectations  of  the  form 
E(X  \d)  can  be  computed  with  high  accuracy.  Let  us  define, 

E  (X  I  d,k)  :  The  conditional  expectation  of  the  random  variable  X,  given  that  the  arrival 
interval  contains  k  packets,  and  has  length  d.  Then, 

OO 

E(X\d)='£tE(X\d,k)e  (B.2.27) 

The  quantities  E(X  I d,k)  depend  only  on  k.  In  Section  B.l  we  show  that  the  quantities 
Lk  =  E(l  I  d,k)  can  be  computed  recursively: 

L0  =  L  i  =  l  (B.2.28) 

Ai2)  aP 

L>=T^L‘-'+-rw;k>-2  (B'2-29) 

where  A„\  i=  1,2,3  can  be  computed  recursively  as  follows: 
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A0)=[l_2-»r1  2“"  jl/xV^jj,  «^3 

^n)=[l-2~"]_1  2-"jn  +  V  A‘2)[||*  n~3 

A^=[l-2~"}|  l+2_n(n+l)+2_"  "]£  A<3)[jj,  n> 3 
A^>  = 1/3,  Ap>  =  2/3,  A$?>  =  7/3 

The  quantities  Z*  =  £  (z  I  d,A:),  can  be  also  computed  recursively: 


(B.2.30) 

(B.2.31) 

(B.2.32) 

(B.2.33) 


Z0  =  0,  Z  i  =  1  (B.2.34) 

a!5)  a!6) 

z‘  =  TT^rz-  +  T^>;t22  (B'2-35) 

and  the  quantities  i  =  4,5,6  can  be  computed  recursively  as  follows: 

A^4)  ;  n  >  3  satisfies  the  recursion  in  (B.2.30) 

A<5)  ;  n>  3  satisfies  the  recursion  in  (B.2.31) 

A<,6)  =  [1  -  2"TVl  +  n2-n  +  n22-n  +  2~""£[ A<6>;  ,  n  >  3  (B.2.36) 

1=2  ^  ^ 

M4)  =  y,  Ai5)  =  y.  46)  =  -y  (B.2.37) 

The  quantities  K*  =  £(/  Id, A:)  can  be  computed  as  follows: 


K0  =  0,  Kj  =  1  (B.2.38) 

Am  A<?) 

rt=7^n-'  +  ^W-k22  <B-2'39) 

and  the  quantities  A^  ,  i  =  7,  8,  9  can  be  computed  recursively  as  follows: 

A(„7)  ;  n  >  3  satisfies  the  recursion  (B.2.30) 
aJ,8)  ;  n  >  3  satisfies  the  recursion  (B.2.31) 

A<’>  =  [1  -2-T‘ (2Z.,-1  +  2-(l  +2LJ  +  n2-(l  +  L._,)  +  7T *£  f \  A?>}  n  >3  (B.2.40) 

1=2  ^  J 


) 
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A?>  =  i,  =  16  (B.2.41) 

From  the  above  formulas  we  see  that  a  finite  number,  M,  of  terms  from  the  infinite 
series  in  (B.2.27),  can  be  computed.  Also,  for  large  k  values,  and  based  on  the  recursive 
expressions,  simple  upper  and  lower  bounds  on  E(X  I  d,k)  can  be  developed.  Those 
bounds  can  be  used  to  tighdy  bound  the  sum 

£  E{X\d,k)e~u^f- 
1  k'- 

Remark  It  can  be  also  proved  that  E  (y  I  d)  =  — — . 

B.3  Interdeparture  Distribution  Analysis 

We  first  give  some  definitions. 

4,m:  Given  k  packets  with  counter  values  equal  to  1  and  m  packets  with 

counter  values  equal  to  2,  the  number  of  slots  needed  by  the  algorithm 
until  the  end  of  the  first  successful  transmission,  after  the  k-multiplicity 
collision  has  been  observed. 

nk,s’-  The  number  of  length  s  interdeparture  intervals  within  a  CRI  wich  starts 

with  a  k-multiplicity  collision.  If  the  length  from  the  first  slot  of  the  CRI 
until  the  first  successful  transmission  is  equal  to  s,  then  this  interval  is 
included  in  the  counting. 

hj:  Starting  with  a  CRP  at  which  the  lag  equals  d,  d>  1,  and  which  follows  a 

successful  transmission,  the  number  of  slots  needed  by  the  algorithm  to 
reach  the  first  CRP  at  which  the  lag  is  one,  and  which  follows  a  slot  con¬ 
taining  a  single  transmission. 
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m^y.  Starting  with  a  CRP  at  which  the  lag  equals  d,  d>  1,  and  which  follows  a 

successful  transmission,  the  number  of  length  s  interdeparture  intervals 
until  the  first  CRP  at  which  the  lag  is  one,  and  which  follows  a  slot  con¬ 
taining  a  single  transmission.  If  the  length  from  the  first  slot  of  the  CRI 
until  the  first  successful  transmission  is  equal  to  s,  then  this  interval  is 
included  in  the  counting. 

P(k,l,h\d):  Given  an  arrival  interval  of  length  d,  the  probability  that  there  are  k 
arrivals  in  it,  that  4,  o=5,  and  that  it  takes  /  slots  for  its  resolution. 

Pk(l)'-  Given  a  k-multiplicity  collision,  the  probability  that  it  takes 

/  slots  for  its  resolution. 

We  also  have: 

H  =  E{h\}  ,  C  =  \H 
c, 

ln(s)}-E{mhs } 

n= 1 

Recursions 

From  the  operation  of  the  algorithm  we  have  the  following  recursions: 

h,m  =0\Vm  ,  P(lkm=0)=0  ;  V  lc>2,  Vm 

fl  ,  ifm=l 

/<),/n  =  l+4i,0  •  Pd0,m  =  1)=|0,  ifm*l 

4 ,m  =  i  +4 ;  WP  \  2~k,  k>2 


(0 


I 
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P(km=S)~ 


1  ,  if  k=l  and  ^=0 
P  On,  0  =5-0  .  if  ^=0,  5>1 


2  1  I  ,1  ^  Vi.m+k-i=S-l),  lfk>2,  s>  1 
1=0  ^  J 


_  Jl  ,  if  5  =  1 
n'-s=  I  0,  if  5*1 


fnjk- is  ,  w.p.  P(lK0*s-l) 
k>2  ;  nu=  j  i  +nk_u  ,  w.p.  />(/*,  0  =5- 1) 


I  /’(/i,  0=5-1)  ,  if  5  >  1 

A \stE  (nksi=  '  k 

l  +  S/>a,.o=0)=l,  if 5=1 

i=2 


*=°’ 1  :  ^  j i  orLLe  •  ^(/)  =  ‘ P(/2.o  = /-2)’  for 


fc>2  ]  i-k- 1 

/>*+2f;/>*(/)=  Z  P«k.o=s)Pk-x  (/-5-1) 


Therefore  for  X  denoting  the  Poisson  intensity, 


/>(Upl</)  =  «-W^-P(/u  =  p)/,H  (/-P-1) 


I 
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Recursions  for  h d 


The  operation  of  the  algorithm  yields  the  following  relationships: 


1  ,  w.p.  Xde  x 
d<A  ;  hd  =  l+/ii  ,  w.p.  e~ 


l+h, ,  w.p.  Ze-™  Pk(l),  l> 2 
*  = 2  K‘ 


d> A  ;  hd  =  l+hd.A+l  ,  w.p.  £e~u  Pk(l) 

*  = 0  *• 


After  taking  expected  values  we  obtain: 


HdiE{hd}=  e~u  +  E{l\d}+e~uHx  +  £  Pk{l)Hr,  d< A 

*>2  l >2  k- 

Hd  =  E{1 1 AJ  +  £  ^  Pk(DHd.A+l  ;  d  >  A 

*20/21 


where. 


=  z  z 

*20/21 


Recursions  for  m„ 


The  operation  of  the  algorithm  yields  the  following  relationships,  where  [J  denotes 


integer  part,  : 


I 
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r 

n i,s ;  w.p.  e  -u\d 

"u  +  :  WP  e~u^j^-Pk(.l)  ;  k>2 


d<A\  m = 


I 

nn-u+m^  .w.p.  P(0ld)  n> 0  /P"(0I  l)/»(/fc,/,p !  1) ;  *>2 
n  +p-t-2« 


1  +  rt*_u  +  m,.,  ;  w.p. 

1  ;  wp.  \e-Md*s-\)  .  f>2 


I 

n>0  «-x*/,(/to=P)/Vi(/-p-l);*22 
n+p+2=i 


For  d  >  A: 


"»<t*  =  "t,  +  ^-a+u  ;  w.p.e  u  />*(/)  ;  *>1 


(W  ry  fl  _\  n  ✓»  ^  i\  t-vi  :r  I  ^  ^ 


=  «*-u  +  fnd-n(A-]y*i.s  ;  *■/>•  X  /5at,o=P) /\-,(/-p-l).  *21.  if  L  r^J  >  1 


_ _ I 


=  I-*-  +  ;w.p.  £  =  (/-r- +«).4>1, 


:,/L  A=fJ  S1 


nk-\.s  +  m,j  ;  W.p.  <? 


d-A 
A- 1  J 


XW_l  A=fJ(A_1))l 


*! 


I  P (ko  =  p)Pk-\0-p-l)  \  k>2 


d- A  . 
A-l  J 


=  1  +/!*-!,,  +m(j  ;  w.p.  e 


*  f 

X 

A-l 

. 

‘H-£if-|(A-1)l 


k ! 


A-l 


■,k>  2 


«*-u  +  «u ;  wPe 


a'1  n> 0 


"+2+p4{-^-Y-J 


e  JTP</*.o  =  P)/W/-p-l);*22 


=  1  +"*-ij  +  "«/.,  ;  w.p.e 


<*-A, 
A-l  ' 


Ze  vp 

nZO  K- 


*-»  A  dI  i  .  .  ..  i  i  ^  4i 

<t.o  -2-t  — — J 


A-l' 
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I 
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+  e  J  A^X(d+p-pA),i  P  }J.d+p-pA).s-\  S  ^l,s  P\(d+p-pA)d) 

l  ‘23 


,-Ud+p) 


l-e 


-x 


/2i3 


i/<p-l>(  "  —  ig^  Uu.WWl]  +  I  I 

1  1~£  L  J  l<nSp  /2l 


Md 

-n(A-iyt/,s  Pte(l) 


+  t/(p-l)£/(s-2)e 


-  U(s-2)Xe~Md+p) 


_0V-X(5-1) 


s-2 

s 

m=j-l-min(p,  $-1) 


,VA  e 


X/7i  n 

^  ^XA.  m 


-X 


(d+p-pA)^  +  , 

\-e~K 


+  U(s-l-p)e  ^  P \(d+p~PA),s-\-p 


+  U{s-2-p)e-*d+s-1)S  e*"  PKm 

m= 0 


Bounds 


For  the  numbers  /Vt  i,  we  used  the  following  bounds: 


0<Nu<lc-\  ;  Vs 
For  the  numbers  Hd,  we  can  prove  that. 


f 


where. 


ot/d  +  P/  ^  Hd  <  CLU  d  +  PB  ,  d  S  1 


I 


oc^a^A-fY/lAyr1  Ff/iA; 


P/=l^d<AeW)  ’  P«  =  max 


-au 


where: 


I 


» 
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I 


Q(d)  = 


E{l\d}  +  au 


E{l\d}-d-Xde 


The  following  simple  bound  on  Mds  has  been  used  : 

0 <Md'S  <H%  =  aud  +  $u 
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APPENDIX  C 


Here,  we  provide  recursive  formulas  for  the  computation  of  Clearly, 


.q  ,q  -|0  otherwise 

Let  /„  be  the  number  of  slots  needed  for  the  resolution  of  multiplicity  n>0  conflict.  Then, 
from  the  operation  of  the  algorithm  we  conclude  that 


p(m)  - 
r  1  ~ 


p(n  + 1) 
rq-\ 


1  fm-fl  1 

W.D.  —  - - 

2  1  n  J  2m_1 


The  upper  pan  of  (C.2)  is  derived  by  considering  the  event  that  the  packet  under  con¬ 
sideration,  together  with  n  of  the  rest  m-1  packets  retransmit  immediately  after  the  initial 
collision.  The  lower  pan  of  (C.2)  is  derived  by  considering  the  event  that  the  packet 
under  consideration  does  not  transmit  immediately,  while  the  n  of  the  rest  m-1  packets 
retransmit  immediately  after  the  initial  collision  and  it  takes  /  number  of  slots  to  resolve  a 
collision  of  multiplicity  n.  The  probabilities  p(ln=l)  can  be  computed  by  similar  reason¬ 
ing.  Averaging  in  (C.2),  we  finally  have  the  following  recursive  formulas  for  Pqm)  for 


P(qm)  =  2-m  "z  P(q-Mmn  ]  +  2-m  l’  Z  \mn  1  P«n=D  (C.3) 

n=0  L  J  n=0  l=\  L  J 
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APPENDIX  D 

We  introduce  a  simpler  notation.  Given  fkj}\<j<3,  we  define: 

Xi^(Xi+X3)A,  i=l,2 

^(X3A)I2  (D.l) 

M-!  —  ^-i+^-3  i  1=1,2 

i=l,2,  ;  P,(e,m):  The  probability  that  T^+e  is  the  m-ih  colli¬ 

sion  resolution  point  after  T)H\  for  the  RAA  in 
cluster  i.  Note  that  m=l  refers  to  the  end  of 
the  first  CRI,  which  starts  at 

i=l,2  ;  P{t)(e,m):  The  probability  that  T^l j  -  T^  =  e,  and  there 

are  m  CRIs  for  the  algorithm  in  cluster  i.  in 

[7lf\  T'ri+l  ]■ 

i=l,2  ;  P{l)(m):  The  probab’lity  that  there  are  m  CRIs  for  the 

algorithm  in  cluster  i,  in  (T(/\  T^h  ]■ 

lx:  The  length  of  a  CRI.  when  it  resolves  all  the 

packets  that  have  arrived  during  an  interval  of 
length  d.  The  number  of  arrivals  per  slot  is 
Poisson  distributed  with  parameter  A..  and 
hd=x. 

0 <n<k  ;  The  number  of  slots  needed  by  either  one  of 

the  two  RAAs  in  the  system  to  transmit  k 
packets,  when  n  of  them  have  counter  values 
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equal  to  one,  and  the  remaining  k-n  packets 
have  counter  values  equal  to  two. 

0<l<e-ll  \  Pi{e,m;l,n): 

,=1-2  J  The  probability  that:  Tjf'+e  is  the  m-th  colli¬ 

sion  resolution  point  after  T%\  for  the  RAA  in 
cluster  i,  T^i  -T(*)>e,  and  the  last  before 
T^+e  collision  resolution  point  for  the  RAA 
in  the  other  cluster  occurs  at  T^+  /,  and  is  the 
n-th  such  point  after  Observe  that 

7t>  +  e  cannot  be  a  collision  resolution  point 
for  both  RAAs  since  we  then  should  have 
T„l\  =r<s)+e. 

t,  =  1,2  ;  Pt jc(e,m  \e,n): 

tc=  1,  if /=2  The  probability  that:  T^l \  =T{n)+e,  T„l\  is 

2,  if /=1 

the  m-th  collision  resolution  point  after 
for  the  algorithm  in  cluster  i,  and  it  is  the  n-th 
such  point  after  T(„s),  for  the  algorithm  in  clus¬ 
ter  ic. 

From  the  operation  of  the  adopted  RAA  we  conclude  the  following: 


P(h.k-\  =m)  =  P(lk- t.o  =m- 1) 

P(lo.k=m)  =  P(lk.o  =  m- 1) 

P({\.  o  =  1)  =  PU  0.0  =  1)  =  Pdo.i  =  2)  =  /»(/,.,  =  2)  =  1 

l^+l}  :  PVik-l=m)  =  2-2 ( P(ko=m-2>+P(lzk-2=m-\)+2P(lk-uo)  =  m- 2) 


Ill 


k>n 
n>  3 

m>2k-l 


;  P(ln.k-n=™)= 2“ 


^(4,  0  =m-2)+P0n.k-n=m -\)+nPUk- l.0~m~2)+ 


CD-2) 


Taking  expectations  with  respect  to  the  number  of  packets  in  the  examined  interval 


we  obtain,  where  |_  j  denotes  integer  part,: 


P(lx=m)=  £  e 
0<k<[  2»±Lj 


P(lk0=m) 


(D.3) 


Furthermore,  it  can  be  proved  that  the  following  recursive  expressions  are  true: 


i=l,2 

e>m 


;  P ,(e,m)  =•{ 


P{lr-y=e)  ;  if w  =  l 


e-\ 


X  Pi(k,m-l)POx,  =e-k)\  if m>2 

i=m- 1 


For, 


A)1  ’if/=2  ,  A  f  if -tS0 


'c  =  ]2,  if /=l 


U(x) 


=  10,  if -r<0 


we  can  also  prove  the  following: 


(D.4) 


(D.5) 


1  =  1,2,  m>  1,  e>m-,  Pi(e,m;0,0)  =  Pi(e,m)P(lXiry>e)  (D.6) 

i=\,2,m>l,k>\,eim\Pl(e+k,  l;e,m)=Plr(e,m;0,0)P(lx^ >k)P(l^_y^e+k)P''1  (iJtv_>>e)(D.7) 


m>2,  e>/7i( 

t=l,2  |  'PtiSe’m:e,  1) =Pl(e,m)P(l^y=e) 


(D.8) 


m>2 

m  i 

e>m 

ten 

1=1,2 


f  :  F.fe.tti ;/,«)=  J)  Pl(e-k,m-l\l,n)P(lx=k)P{lx^>e-f)P'\l^  >e-k-h  + 

1  *=i 


+U(l-m)‘ ^  P,t{Un\e-k,m-\)P{.llt=k)P{lK>e-[)p-\lx>>l+k-e) 

k=e-i+l 


(D.9) 
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e  -1  min  (e  -7— 1 ,  *  1  ~ m ) 

■,P,t<(e,m\e,ny=  £  £ 


/=«+!  *=1 


./>(/,  =e-t)P~\lx  >e-k- /)  + 

\  «e 


t-1  min(«-/-l,«+l-/i) 

+  I  £  P.'ie-^n-lMm-nPd^PV^-DP-'U^e-k-O  (D.10) 

/=m-l  *=1 

We  therefore  obtain: 


p(‘\e,m)  =  £  Pij'  ( e,m\e,n ) ;  t'=l,2 

/I 


(D  .11) 


P(i\m)  -  £  />(,)(e,m) ;  t=1.2 

<f  >m 


(D.12) 


The  mean  of  the  distribution  in  (D.  12)  is  the  quantity  E{Ni  ( (Kj)  ,  A ) }  in  (V.3). 
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